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1 Introduction 



Present measurements of the vector boson-fermion couplings at LEP and SLC accurately con- 
firm the Standard Model (SM) predictions at the 0.1 - 1% level which may readily be 
considered to be evidence for the gauge boson nature of the W and the Z. Nevertheless the 
most crucial consequence of the SU{2) x U{1) gauge theory, namely the specific form of the 
non-Abelian self-couplings of the W, Z and photon, remains poorly measured to date. A direct 
and more accurate measurement of the trilinear self-couplings is possible via pair production 
of electroweak bosons in present and future collider experiments (py^M^" at LEP2, W'y, WZ 
and W^W~ at hadron colliders). 

The major goal of such experiments at LEP2 will be to corroborate the SM predictions. If 
sufficient accuracy is reached, such measurements can be used to probe New Physics (NP) in 
the bosonic sector. This possibility raises a number of other questions. What are the expected 
sizes of such effects in definite models of NP? What type of specifically bosonic NP contribu- 
tions could have escaped detection in other experiments, e.g. at LEPl? Are there significant 
constraints from low-energy measurements? Although we shall address these questions, the 
aim of this report is mostly to elaborate on a detailed phenomenological strategy for the direct 
measurement of the self-couplings at LEP2, which should allow their determination from data 
with the greatest possible accuracy. 



2 Parametrization, Models and Present Bounds on TGC 

We shall restrict ourselves to Triple Gauge boson Couplings (TGC) in most of the report (possi- 
bilities to test quartic couplings at LEP2 are extremely limited). Analogous to the introduction 
of arbitrary vector and axial-vector couplings gy and qa of the gauge bosons to fermions, the 
measurements of the TGC can be made quantitative by introducing a more general WWV 
vertex. We thus start with a parametrization in terms of a purely phenomenological effective 
Lagrangian Q 0, ^ [V = •y or Z] 



^'-eff — 9WWV 



gXV^ {W-W+'' - W+W-") + Ky W^W~V^' + (1) 



which gives the most general Lorentz invariant WWV vertex observable in processes where the 
vector bosons couple to effectively massless fermions. Here the overall couplings are defined as 
gww-y = e and gwwz = e cot 6w, W^y = d^JWy - dyW^, and V^y = d^V^ - dyV^,. For on-shell 
iWe use £0123 ^ i. 
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photons, Qiiq'^ = 0) = 1 and gjiq"^ = 0) = are fixed by electromagnetic gauge invariance f\ 
Within the SM, at tree level, the couplings are given hj gf = gj = = t^-y = ^, with all other 
couplings in (0) vanishing. Terms with higher derivatives in (|I|) are equivalent to a dependence 
of the couplings on the vector boson momenta and thus merely lead to a form-factor behaviour 
of them. We also note that g^ , and Ay conserve C and P separately, while g^ violates C 
and P but conserves CP. Finally g^ , ky and Xy parameterize a possible CP violation in the 
bosonic sector, which will not be much studied in this report, as it may be considered a more 
remote possibility for LEP2 studies 0. However, there exist definite and simple means to test 
for such CP violation, see section The C and P conserving terms in C^jf^ correspond to the 
lowest order terms in a multipole expansion of the VT— photon interactions: the charge Qwi 
the magnetic dipole moment nw and the electric quadrupole moment qw of the 0: 

6 6 

Qw = egj , nw = {gj + K,^ + X-y) , qw = 2~ (^7 " -^7) • (2) 

For practical purposes it is convenient to introduce deviations from the (tree-level) SM as 

Agff = ((yff — 1) = tanOw^z , = (k^ — 1) = , (3) 

Akz = {kz - 1) =tan9w{xz + 6z) , 
Xy =yy, Xz =tan6wyz- 

For completeness (and easy comparison) the correspondence of the most studied C and P 
conserving parameters has also been given for another equivalent set {6z, xy, yv), which was 
used in some recent analyses ||^, 0] . 



2.1 Gauge-invariant Parametrization of TGC 

Any of the interaction terms in (|I|) can be rendered SU{2) x U{1) gauge invariant by adding 
to it interactions involving additional gauge bosons [||, and/or additional Would Be Goldstone 
Bosons (WBGBs) and the physical Higgs (if it exists) @, |ll|] . However, one needs to consider 
SU{2) X f/(l) gauge invariant operators of high dimension in order to reproduce all couplings 
in (|l|). For example, if the Higgs particle exists one needs to consider operators of dimension 
up to d = 12. Depending on the NP dynamics, such operators could be generated at the NP 
mass scale A^p, with a strength which is generally suppressed by factors like {mw / A^pY''^ or 
{y/s/ AisfpY~^ |12|, Accordingly, the gauge invariance requirement alone does not provide 
any constraint on the form of possible interactions. Rather it is a low energy approximation, 
the neglect of operators of dimension greater than 4 or 6, which leads to relations among the 
various TGCs. 

Such relations among TGCs are highly desirable, given the somewhat limited statistics 



accessible at LEP2. They were first derived in [|1^, § by imposing approximate global SU{2) 



^For ^ deviations due to form factor effects are always possible, see section p!4| below in this connection 
^Data on the neutron ele( 
at the 10"'^ level is accepted 



•^Data on the neutron electric dipole moment allow observable effects of e.g. at LEP2 only if fine tuning 
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symmetry conditions on the phenomenological Lagrangian ([1|). In the next subsection we 
present them following an approach based on SU{2) x U{1) gauge invariance and dimensional 
considerations. The connection to the approach based on "global SU{2)" symmetry will be 
discussed at the end. 

In order to write down all allowed operators of a given dimensionality one must first identify 
the low energy degrees of freedom participating in NP. We assume that these include only 
the SU{2) X U{1) gauge fields and the remnants of the spontaneous breaking of the gauge 
symmetry, the WBGBs that exist already in the standard model. If a relatively light Higgs 
boson is assumed to exist, then NP is described in terms of a direct extension of the ordinary 
SM formalism; i.e. using a /inear realization of the symmetry. On the other hand, if the Higgs 
is absent from the spectrum (or, equivalently for our purpose, if it is sufficiently heavy), then 
the effective Lagrangian should be expressed using a nonlinear realization of the symmetry. 



2.1.1 Linear Realization 



In addition to a Higgs doublet field $, the building blocks of the gauge-invariant operators are 
the covariant derivatives of the Higgs field, -D^$, and the non-Abelian field strength tensors 
W^/^i/ = — gW^ X Wi, and -B^^, of the SU{2)l and U{1)y gauge fields respectively. 

Considering CP-conserving interactions of dimension d = 6, 11 independent operators can 
be constructed [|r^, ^ |TD[. Four of these operators affect the gauge boson propagators at tree 



level [|l^] and as a result their coefficients are severely constrained by present low energy data 



IP, |TOl- Another subset of these operators generates anomalous Higgs couplings and will be 



discussed in section |10.4| below. Here we consider the three remaining operators which do 
not affect the gauge boson propagators at tree-level, but give rise to deviations in the C and 
P-conserving TGC. Denoting the corresponding couplings as aw<p, 0:34,, and aw, the TGC 
inducing effective Lagrangian is written as 

(4) 

with g, g' the SU{2)l and U{1)y couplings respectively. Replacing the Higgs doublet field by 
its vacuum expectation value, $^ (0,f/-\/2), yields nonvanishing anomalous TGCs in (Q), 

2 

^91 = , = -^{Akz - Agf) = aw4, + as,), , = \z = aw , (5) 

Cw -^w 

where sw = sin^^, cw = cos 6*1^. The normalization of the dimension 6 operators in has 
been chosen such that the coefficients ai correspond directly to Ak^ and A^. It should be noted 
that, as the NP scale Aatp is increased, the ai are expected to decrease as {mwf-^Np)'^ ■ 

This scaling behaviour can be quantified to some extent by invoking (tree-level) unitarity 
constraints [|17|, |18|, |T9|, |T^ . A constant anomalous TGC leads to a rapid growth of vector boson 



pair production cross-sections with energy, saturating the unitarity limit at ^/s = Ajj. A larger 
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value of Kjj implies a smaller TGC For each of them the unitarity relation may be written 
as 



17, m 



,2 / \ 2 

my/ \ I I / Tnyv 

\aw\ ^ 19 1 -; — I , \aw4,\ ^ 15.5 



Ar 



Ar 
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XT 



(6) 



For any given value of the corresponding scale Kjj provides an upper bound on the NP scale 
Atvp- Conversely, a sensitivity to small values of an anomalous coupling constant is equivalent 
to a sensitivity to potentially high values of the corresponding NP scale. Applying for 
A[/ = 1 TeV, we get \aw\ — 0.12, \aw4>\ — 0.1, |aB</,| — 0.3. Since these values are larger 
than the expected LEP2 sensitivity by less than a factor 3, it is clear that LEP2 is sensitive to 
Aatp ^ 1 TeV. Thus a caveat is in order: for these low values of A^rp the neglect of dimension 
8 operators may no longer be justified, leading to deviations from the relations ® [10 • 



2.1.2 Nonlinear Realization 

In the absence of a light Higgs a non-linear approach should be used to render -C^]^^ gauge 
invariant. The SM Lagrangian, deprived of the Higgs field, violates unitarity at a scale of 
roughlyQ47rf ~ 3 TeV, so that the new physics should appear at a scale Aatp < Attv. Technically 
the construction of gauge-invariant operators follows closely the linear case above, except that 
in place of the scalar doublet $ a (unitary, dimensionless) matrix U = exp(icJ -t/v)^ where the 
uji are the WBGBs, and the appropriate matrix form of the SU (2)^ x f/(l)y covariant derivative 
are used. The so-called "naive dimensional analysis" (NDA) [l2| dictates that the expected 
order of magnitude of a specific operator involving h WBGB fields, d derivatives and w gauge 
fields is ~ f ^A^P {l/vf {1/KnpY {g/^NpT- Applying NDA to the terms in Eqs. (1), we see 
that and Akv are of 0(m^/A^p). In other words, just as in the linear realization, these 
terms are effectively of dimension 6 (in the sense that there is an explicit factor of 1/A^p). 
On the other hand, we see that the W^^W'^uV'' term is effectively of dimension 8, i.e. the 
coefficient Xy is expected to be of order m^/A^p. Thus, within the nonlinear realization 
scenario, the Xy terms are expected to be negligible compared to those proportional to AgY 
and Akv- Accordingly there remain three parameters at lowest dimensionality, which can be 
taken as g^ , kz and k^. 



2.1.3 Operators of Higher Dimension and Global Symmetry Arguments 



As mentioned in section |2.1| , one may argue that relations like in would not even be approx- 
imately correct if A^vp is substantially smaller than 1 TeV, since higher dimensional operators 
are no longer suppressed, and may even be more important than the dim = 6 operators |20 . 



In fact, as far as the 5 C and P conserving TGC in (|l]) are concerned, the most general choice 
can be realized by invoking two dim = 8 operators in addition to the 3 terms in (^) [|lO, |ll], |2^ . 

One should caution that this estimate of Ajvp follows directly from analogy with low energy QCD and 
Chiral perturbation theory |2l|] , where v = f^^ and A ~ Alp are known, while in the present context Ajvp is 
essentially unknown. It should be taken as a rough order of magnitude estimate only. 
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Requiring restoration of an SU{2) global ("custodial") symmetry for g' ^ (i.e in the limit of 
decoupling B field) implies |Q the coefficient of one of these two operators to vanish, because 
it violates SU (2) globalQ independently of the B field. In that way one recovers the constraints 
between Ak^ and Akz in in both the nonlinear realization and in the linear realization 
at the dim = 8 level. Nevertheless a second dim = 8 operator spoils the relation, = Xz in 
d^). One may neglect this term (which vanishes in the limit g' — > 0) by imposing exact SU{2) 
at the scale A^rp, which in our context is similar to neglecting the tt^ — vr" mass difference in 
strong interaction physics. 

Largely these are simplifying assumptions only, intended to reduce the number of free pa- 
rameters. Motivated by the previous discussion we recommend two sets of three parameters 
each for full correlation studies between anomalous couplings at LEP2: 



setl = {Agf, Ak^, Akz) with = Xz = 0. These correspond to the operators of 
lowest dimensionality in the nonlinear realization. A reduction to 2 parameters (using 

2 

Ak^ = —^{Akz — Agf)) is achieved by assuming p, ^ custodial SU{2) for g' ^ 

set2 = (Agff, Ak^, A^) with Xz and Ak,z given by (|^). It is this set which has been 
used in this report for the determination of precisions achievable from WW production 
at LEP2, presented in sections 5-9 as limits on the parameters a^^, cxw^ and aw defined 

by m. 



Expressing results in terms of Agf , An^, etc. will be useful for ease of comparison with 
published hadron collider data p5|, pG] , 

In addition, it would clearly be of interest to present fits to each of the parameters in 
in order to reduce the dependence of the analysis on specific models. However, this can only 
be achieved bearing in mind the limited data which will be available from LEP2, and the 
correlations inherent in the extraction of many parameters from the data. We return to this 
point in sections 3.1, 4.2 and 5.1 below. 



2.2 Present constraints on TGC 



The errors of present direct measurements, via pair production of electroweak bosons at the 
Tevatron, are still fairly large. The latest, best published 95% CL bounds by CDF and DO are 



obtained from studies of W'y events ^ 



-1.6 < Ak^ < 1.8 , -0.6 < A^ < 0.6 

but constraints from the study of WW, WZ — >• iujj, ££jj events are becoming competitive 
and should lead to 95% CL bounds of roughly —0.65 < Ak^ < 0.75, |A^| = \Xz\ < 0.4, once the 

^Note that there is no contradiction with the SU{2)l x U{1)y local invariance of all these operators, since 
SU{2) custodial is a different symmetry from the SU{2)l global ||2^ . 
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already collected run lb data are fully analyzed. Increasing the integrated luminosity to 1 fb ^ 
with the Fermilab main injector is expected to improve these limits by another factor 2 . 



Note that these latter bounds assume the relations between anomalous couplings as given by 
(H) with aw(j, =OiB(i,- In addition, the Tevatron measures these parameters at considerably 
larger momentum transfers than LEP2 and, hence, form factor effects could result in different 
measured values at the two machines. 

Alternatively, constraints may be derived also from evaluating virtual contributions of 
TGC to precisely measured quantities such as {g — 2)^ the b S'j decay rate [^, |3( 



B K^*^fi~^fi~ the Z bb |^ rate and obhque corrections]^, |T0[ {i.e. corrections to 
the W, Z, 7 propagators). Oblique corrections combine information from the recent LEP/SLD 
data, neutrino scattering experiments, atomic parity violation, //-decay, and the W^-mass mea- 
surement at hadron colliders. 

When trying to derive TGC bounds from their virtual contributions one must make assump- 
tions about other NP contributions to the observable in question. In the linear realization, for 
example, Higgs contributions to the oblique parameters tend to cancel the TGC contributions 
and as a result the TGC bounds are relatively weak for a light Higgs boson [|l^. In general, 
there are other higher dimensional operators which contribute directly to the observable, in 
addition to the virtual TGC effects. Bounds on the TGC then require to either specify the 
underlying model of NP completely or to assume that no significant cancellation occurs. The 
bounds on the TGC parameters in (|^) due to virtual effets thus depend on the underlying 
hypotheses and are of 0(0. 1) to C(l) g |TU|, 

More stringent bounds are obtained by comparing the higher dimensional operators which 
induce TGC with those operators which directly induce oblique effects (see Section 10.1). In 
simple models the coefficients of these two sets of operators are of similar size and hence the 



stringent LEPl bounds on the latter ||3J] indicate that one should not expect anomalous TGC 
above (9(0.01). One should stress, however, that no rigourous relation between oblique effects 
and TGC can be derived except by going to specific models of NP. Therefore, these stringent 
bounds must be verified, by a direct measurement of the TGC at LEP2. 



2.3 Virtual Contributions to TGC in the MSSM Q 

Definite TGC contributions are certainly present at the loop level in any renormalizable model, 
although such loop effects contribute to TGC with a factor of ((/^/IGtt^) ^ 2.7 10~^, being 
therefore too small a priori to be observed at LEP2. For instance, SM one-loop TGC predictions 



are known [||, |3|, ^ and give, at = 190 GeV, Ak^ (Akz) ^ 4.1-5.7 lO^^ (3.3-3.1 lO^^), 
for rriHiggs = 0.065-1 TeV and m^p = 175 GeV |3^. (Contributions to Xy are about a factor of 
3 smaller). One may, however, expect that the "natural scale" {g'^/lGn'^) could be substantially 
enhanced if, for example, some particles in the loop have strong coupling and/or are close to 

complementary study of virtual MSSM contributions to the e^e" VF+VF^cross-section is done in the 
New Particle chapter of these proceedings. 
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their production threshold. To obtain a "reference point" it is thus important to explore more 
quantitatively how far one is from the LEP2 accuracy limit, within some well-defined model of 
NP. We here use the contributions of the (MSSM) as an example. These contributions were 
calculated independently by two groups in the framework of the Workshop. We summarize the 
main results, referring for more details to refs. 13^, 



SUGRA-GUT MSSM 


Unconstrained MSSM (maximal effects) 


Ao, mo, Ml/2 = 300,300,80 (GeV), 
tan/? = 2 (yu<0) 


tan/3 =1.5; M^^ ~ 95, 130; Mf ~ 20-132 (GeV) ; 
mH+ ^ 95; mj>, ~ 45; m,- ^ 92-110; rriq ~45-800 GeV 


|Ak^| = 0.44 10-2, \ Akz\ = 0.72 IQ-^ 


Ak^ = 1.75 10-2, Akz = 0.84 lO'^ 



Table 1: Ak^^z (as defined in eq. |I|) in MSSM at ^/i = 190 GeV. (Contributions to Ay are 
about a factor of 2-3 smaller). 



Naively, TGC are obtained by summing all MSSM contributions to the appropriate parts in 
eq. (|I|) from vertex loops with entering 7 (or Z) and outgoing VT"*", W~ . But as is well-known, 
the vertex graphs with virtual gauge bosons need to be combined with parts of box graphs 
for the full process, e^e~ — W^W~, to form a gauge-invariant contribution. The resulting 
combinations define purely s-dependent [] TGC p6[. In table |l] we illustrate our results for 
(s-dependent) contributions in two different cases. First, for a representative choice of the free 
parameters in the more constrained MSSM spectrum obtained from the SUGRA-GUT 
scenario 0]: the only parameters are the universal soft terms mo, M1/2, Aq at the GUT scale, 
tan/3 (and the sign of fi). Second, we give one illustrative contribution, obtained from 
a rather systematic search of maximal effects in the unconstrained MSSM parameter space. 
The largest contributions are mostly due to gauginos and/or some of the sleptons and squarks 
being practically at threshold. One may note, however, that some individual contributions, 
potentially larger, were quite substantially reduced when the present constraints on the MSSM 
parameters are taken into account |Q. Even these maximal contributions hardly reach the 
level of the most optimistic accuracy limit expected on TGC (compare section ^ below). One 
should also note that radiatively generated TGC generically have a complicated ^/s form factor 
dependence as well as contributions from boxes, which are well approximated by an expansion 
in 1/Amp only when one probes well below threshold. 



2.4 TGC from extra Z' Q 

A light and weakly coupled Z' provides an illustrative example of relatively large deviations 
of the TGC from their SM values and of strong form-factor effects Consider an extra 

^By definition, t and u-dependent box contributions are left over in this procedure. We have evaluated |]38| a 
definite (gauge-invariant) sample of this remnant part, the slepton box contributions, and found them negligible, 
~ 0.1 ig'^/lGn'^) ~ 3 10"'' at most, at LEP2 energies. 

®A complementary study can be found in the Z' working group chapter of these proceedings 
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gauged U{iy symmetry with associated coupling g[, whose vector boson Z' is relatively light, 
say Mz' — 200GeV. For such a boson to remain undetected at LEPl and CDF, it must have 
rather small couplings to fermions: A = sm{6w)gi/gi < 0.2 or less 0]. However, this new 
Z' might be only part of the new physics beyond the SM, and we parametrize this by gauge 
invariant higher dimensional operators. For illustration, let us focus on the dim = 6 operator 



■ TC 



(7) 



where B'^^^ is the new U'{1) field strength. This operator has a part linear in Wf^ inducing 
unusual mixing through the kinetic terms, from which LEPl data put upper bounds on A and 
e. The other piece is quadratic in and brings anomalous contributions to py-pair production 
at LEP2, which may be enhanced at will by approaching the Z' pole. Within a gauge-invariant 
framework, enlarging the symmetry group has given us enough freedom to escape the more 
stringent LEPl constraints on the coefficient of the similar operator Obw H, HO) of Eq. (^9]) . 
Having such an (admittedly contrived) counter-example to (depending on the 3 parameters 
Mz', A < 0.2 and |e| < 0.2), it is instructive to see how it fits into our TGC parametrization. 

The normal way of extracting the pre- 
dictions of this model for ly-pair produc- 
tion would be to add all the amplitudes for 
e^e^ W~^W~, namely the t-channel u 
pole, and s-channel 7, Z and Z' poles, in- 
cluding the contributions of Ob'w in the lat- 
ter. Alternatively, the correct angular depen- 
dence in e~^e~ — >■ iy+iy~from such a Z' is 
recovered through the introduction of "pro- 
cess - dependent" TGC form factors: the Z' 
exchange only contributes to the J = 1 par- 
tial wave and the TGC of Eq. (|l]) allow to 
parameterize the most general J = 1 am- 
plitude. For the case at hand one can al- 
ways find TGC {gf , , gj, k"') matching the 
Z' parameter dependence in this particular 
ee-WW channel, but these TGC will depend 
on the incoming electron's coupling to the Z 
and the photon. 




^gf 

The deviations Agf vs. AgJ for ^/s = 205GeV 
and Mz' = 210GeV. For each A ranging from 
(plain curve) to 0.2 (smallest dashes), e is 
limited to satisfy today's W mass accuracy, 
\SMw\ < IGOMeV (LEP2's \5Mw\ < 45MeF 
for the thick curves). 



In general, a non-zero AgJ is needed to match the precise t-dependence, but in such a 
process-dependent approach, this does not imply any violation of charge conservation. Finally 
one should note that the Z' described above would also appear in e~e+ — >• qq, at LEP2 

and thus all channels need to be searched for NP effects. 



3 The W Pair Production Process 
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3.1 Phenomenology of On-shell WW Production 



Deviations of the TGC's from their SM, tree level form are most directly observed in vector 
boson pair production. At LEP2 this is the process e~e~^ W~W~^, which, to lowest order, 
proceeds via the Feynman graphs of Fig. |l]. We start by describing the core process, including 
the W decay into fermion anti-fermion pairs in the zero-width approximation, since most of the 
effects of anomalous couplings can already be understood at this level. A full simulation of the 
signal will, of course, need refinements such as finite width effects, the ensuing contributions 
from final state radiation graphs and the inclusion of t-channel vector boson exchange graphs 
for specific final states such as e'vud. The simulation of this full e^e~ —>■ 4 fermions process 
will be discussed later. It is instructive to consider first the individual contributions of s-channel 



e-{a) 





+ 1 



'^/\/\/v 



w-{x) 



w+Cx) 



Figure 1: Feynman graphs for the process e~^e W^W . 

photon and Z exchange and of t-channel neutrino exchange to the various helicity amplitudes 
for the process e~e+ — >• W~W^ 0, 

M{(y,\,\) = M= M., + Mz + . (8) 

Here the e~ and e+ helicities are given by cr/2 and — o"/2, and A and A denote the W~ and 
helicities. Let us define reduced amplitudes ^A by splitting off the leading angular dependence 
in terms of the c?-functions d'^'^ where Jq = 1,2 denotes the lowest angular momentum 
contributing to a given helicity combination. In the cm. frame, with the e~ momentum along 
the z-axis and the W~ transverse momentum pointing along the x-axis, the helicity amplitudes 
are given by^ 

A4(a, \\:e) = V2a ^.,a,a(^) <:a-a(^) • (9) 

For (A, A) = (±,=f), i.e. |A — A| = 2, only t-channel neutrino exchange contributes and the 
incoming electron must be lefthanded. The corresponding amplitudes are given by 

M{-1, A, A = -A; 6) = -V2e' \ A sinO (1 - A cos0)/2 . (10) 

sm 9w 1 + 13^ — 2(3 cos 

s-channel photon and Z exchange is possible only for |A — A| = 0, 1. The corresponding reduced 
amplitudes can be written as 



'As compared to Rcf. [|j a phase factor (—1)'*' is absorbed into the definition of the W'^ polarization vector. 
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Mz 



z_ 

XX 



1 

2 sin^^ 



w. 



m 



2 ' 



Bx-x~ 



1 + 13^- 2/3cos^ 



C 



XX 



Here s denotes the e+e center of mass energy and (3 = a/1 — 4m^/s is the velocity. The 
subamphtudes , B and C are given in Table ^ 



Table 2: Subamphtudes for Jq = 1 helicity combinations of the process e~e+ W~W^, as 
defined in Eq. (|TT]). /3 denotes the W velocity and 7 = \/sj2mw- The abbreviation = 
gX \ ^^v ^ is used. 



++ (^}^ + 272Ay + ^(Ky + Ay-272Ay) 

— + 27'Av - ^(Ky + Ay - 27'V) 

+0 i{n - igX + PgX + ^(^y - V)) 

0- Tl/l' + ^gX + /^^^s'^ - M'^F - V)) 

0+ lifx + ^^^r - + - V)) 

-0 lUX - wX - PaX - U~^v - Ay)) 

00 gX + 





Cx-x 


cr,A — A 


1 


1/7^ 


-asin^/y2 


1 


1/7^ 


-CTsin^/y2 


27 


2(l + /?)/7 


(1 + a cos 9) /2 


27 


2(l + /?)/7 


(1 + acos^) /2 


27 


2(l-/5)/7 


(1 - a cos 0) /2 


27 


2(l-/5)/7 


(1 - a cos 0) /2 


27^ 


2/7^ 


-cr sin ^/v^ 



One of the most striking features of the SM are the gauge theory cancellations between 7, Z 
and neutrino exchange graphs at high energies. Within the SM the only non- vanishing couplings 
in the table are gi = n = 1 and /s = 2 for both the photon and the Z-exchange graphs. As 
a result AJ^^^ = A^^ and the jSA^ terms in Eq. ([TT| ) cancel, except for the difference between 
photon and Z propagators. Similarly, the Bxx term in Aiu and the ^o-.-i term in Aiz cancel 
in the high energy limit for all helicity combinations. While the contributions from individual 
Feynman graphs grow with energy for longitudinally polarized W^s, this unacceptable high 
energy behavior is avoided in the full amplitude due to the cancellations which can be traced 
to the gauge theory relations between fermion-gauge boson vertices and the TGC's. 

LEP2 will operate close to W pair production threshold and these cancellations are not yet 
fully operative. For example, at i/i = 190 GeV one has /? = 0.54, /? s/{s — m\) = 0.70, and 
1/(3 = 1.87 instead of unity. As a result, the linear combinations of couplings which enter in 
Ai^ and Aiz are quite different from their asymptotic forms. In particular the 7^ enhancement 
factors are still small, the (±, ±) and (0, 0) amplitudes are not yet dominated by individual 
couplings, and interference effects between different TGC are very important. 

Table 1^ shows that only seven W~W^ helicity combinations contribute to the Jq = 1 channel 
and the various WWV couplings enter in as many different combinations. This explains why 
exactly seven form factors or coupling constants are needed to parameterize the most general 
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WWV vertex. Since we have both WW Z and WW'~) couphngs at our disposal, the most 
general J = 1 amplitudes Ml = M{a = —1, A, A) and Mr = M{a = +1, A, A) for both left- 
and right-handed incoming electrons can be parameterized. Turning the argument around one 
concludes that all 14 helicity amplitudes need to be measured independently for a complete 
determination of the most general WW'y and WWZ vertex. 

A first step in this direction is the measurement of the angular distribution of produced 
W^s, da /d cos 6. In terms of the reduced amplitudes A^o-,a,a of ® this distribution is given by 



da 



d cos 9 



7X0^(3 

4s 



E [ 

o-=±l 



sin^e 



\M^,++? + \M,,-\'' + \M<r,m 



'1 + cr cos ( 



(1 



(TCOSI 



4 

+ -(1 + cos^ 



sin^^ 



iA^<x,o+r+i-M^,-or 

2 1 



sin^^VK(l + /?^-2/3cos( 



(12) 



Due to the different ci-function factors amplitudes with different values of A — A can be separated 
in principle. In practice, the additional 6'-dependence of the neutrino exchange graphs (the 
terms in Eq. ([TT|)) distorts these angular distributions and leads to contributions from the 
individual W~W^ helicity combinations as shown in Fig. ^. In fact, the interference with the 
i/-exchange graphs can be used to further separate the various s-channel helicity amplitudes. 
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Figure 2: Angular distributions da/d cos 9 for e^e+ — > W~W~^: SM contributions from fixed 
W-W+ hehcities (AA) at = 190 GeV. 
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Due to the V — A structure of the VT-fermion vertices the decay angular distributions of 
the W^s are excellent polarization analyzers and a further separation of the various W^W~ 
helicities can be obtained |^. These decay distributions are most easily given in the rest 
frame of the parent W . Choose the e~e^ W~W^ scattering plane as the x — z plane with 
the 2;-axis along the W~ direction and obtain the rest frames by boosting along the z-axis. 
In the W~ frame we define the momentum of the decay fermion for W~ —>■ /1/2 as 

Pi ~ (-'-5 sin^i COS01, sin6'i sin0i, cos6'i) , (13) 

and, similarly, for fzfii the anti-fermion momentum in the frame is given by 

P4 = (1, sin6'2 COS02, -sin6'2 sin02, -cos6'2) , (14) 

Thus, 6i = Q corresponds to the charged lepton or the down-type (anti)quark being emitted in 
the direction of the parent W^. 

Neglecting any fermion masses, the W~ ^ i^v decay amplitude is given by 

-Md(A) = ^^ 4(^1, </-!), (15) 
where the angular dependence is contained in the functions 

(£_,£o,^+)(^i,0i)= (^-^(l + cos^i)e-^<^\ -sin^i, -L(l - cos^^) e^^^^ . (16) 

An analogous expression is obtained for the decay amplitude. 

The production and decay amplitudes can easily be combined to obtain the five-fold differ- 
ential angular distribution for the process e~e^ —>■ W~W^ —>■ /1/2 73/4, in the narrow W^-width 
approximation [0, ||, 

d'a (e-e+ ^ ^ fj, hh) ^ ^ (^Wjw f f \ PAW ff) 

d cos e d cos Oi d (f)^ d cos 62 d (t)2 1287rsV8vry ^ ^ 

cr,X,X,\',X' 

X Da,A'(^1,0i) I^A,A'(^-^2,02 + 7r) . (17) 

Here the production amplitudes M{<J, A, A) are given in Eq. (|^) and the Z^a.a' are given by 

Dx,x'i9,<P) = hie,<P)tyi9,<P). (18) 

The information contained in the five-fold differential distribution ( [T7| ) can be used to isolate 
different linear combinations of WWV couplings and hence reduce the possibility of cancel- 
lations between them. For example, by isolating W^W~ pairs which are both transversely 
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polarized (and hence give 1 + cos^^j decay distributions) the combinations gY + 27^Ay are de- 
termined which appear in the production amplitudes and Ai (see Table |]). Similarly, 

longitudinal W's produce a characteristic sin^^j decay distribution. The isolation of LT+TL 
and of LL polarizations of the two W^s allows independent measurements of the combinations 
f3=gY + f^v + and gY + 27^ /ty, respectively, and thus the three C- and P-conserving 
anomalous couplings]^ may be isolated. 

Additional information is obtained from the azimuthal angle distributions of the decay 
products. A nontrivial azimuthal angle dependence arises from the interference between he- 
licity amplitudes for different or different W~ polarizations. The large A^+_ and 
amplitudes, which arise solely from neutrino exchange, can thus be put to use: interference 
with these large amplitudes can amplify the effects of anomalous couplings. 

The observation of azimuthal angular dependence and correlations is particularly important 
for the study of CP-violating effects in W~W^ production 0, The methods suggested in 



section ^ below for TGC determination from data can all be used for this purpose, and the 
reader is referred to the literature for details of procedures using density matrix and opti- 



mal observable analyses. Similarly, the study of rescattering effects between the produced 



W pairs, i.e. the presence of nontrivial phases in the production amplitudes, relies on the 
interference with the phase factors introduced by the azimuthal angle dependence of the de- 
cay amplitudes. We do not explicitly discuss these techniques here but rather refer to the 
literature 0, 



WW decay channel 


Decay fraction 


Available angular information 




/ = e: 14% 
/ = /i: 14% 
/ = r: 14% 


cos6' 


{cos9i,(f)i) {cos 9 j, (pj) folded 


JJJJ 


49% 


cos6' 


{C0s6j-^, 0ji)folded (C0s6'j2, ^jjfolded 


iuiu 


9% 


cos 6* 


(cos 6*1, 01 ) (cos 6*2, 02) 

2 solutions 



Table 3: Availability of angular information in different WW final states. The production 
angle is denoted by 9 and {9ij,(f)ij) denote decay angles for W (leptons, jets) respectively, 
(cos 6'j, 0j) folded implies the ambiguity cos 6'j ^ — cos6'j, 0j ^ 0j + vr incurred by the inability 
to distinguish quark from antiquark jets. 

The application of (|T7|) to experimental data must take account of some restrictions in 
the ability to determine the angles involved: in the case of hadronic W decays, and in the 
absence of any quark charge or flavour tagging procedure, the fermion and anti-fermion cannot 
be distinguished; also, in the case where both Ws decay leptonically, a quadratic ambiguity is 

^°Note however that if relations among TGC such as those in eq. (^) are relaxed, it will not be easy to 
distinguish Kj from kz (or Xj from Xz) with unpolarized beams, since these both feed the same helicity 
amplitudes in table ||. 
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encountered. The ambiguities in each of the three WW final states jjiiy, jjjj and iuiu, where 
j represents the jet fragmentation of a quark or antiquark and (/z/) the products of W decay 
into lepton-antilepton, are summarized in table 0. 



3.2 Four-fermion production and non-standard TGC 



Most studies of TGC so far have been made with zero width simulated data and with an analysis 
program based on the same assumptions. This procedure might neglect some important effects, 
however, and the corresponding physics issues will be discussed in this subsection. These are 
the influence of a flnite W-width, of background diagrams, i.e. graphs other than the three 
W-pair diagrams of Fig. |l|, and the influence of radiative corrections (RC) in particular the 
dominant QED initial state radiation (ISR). 

At the moment there are many Monte Carlo (MC) programs for four fermion production, but 



only two of them can at present study the above issues, namely ERAT0[46| and EXCALIBUR^ 
For a detailed description we refer to the WW event generator report, but we make a few 
comments here. Although the programs can study non-standard TGC effects [^, |49[| for all the 
channels of Table ^, we will only consider the jjiiy case in the following. More speciflcally we 
will study e'VeUd or ^'v^ud flnal states. The amplitude for these flnal states consists of 20 and 
10 diagrams, respectively, of which 3 are the W-pair diagrams of Fig. |l]. Since the four fermions 
are assumed to be massless in the calculations, cuts have to be applied to avoid singularities 
in the phase space. Experimental cuts usually have this effect as well. In the case of only 
three diagrams such cuts are not required. ISR is incorporated following the prescription of 



Ref. [|50|. In table 0, we list a number of differential cross-sections which have been calculated. 



Standard Model 


non-standard TGC 


physical assumptions 


(7 SM,on 


0'AN,on 


rw = o 


0'SM,off 
0" SM, off, cuts 


^AN,off 
CAN,off,cuts 


3 diagrams 


C" SM,all 


O' AN, all 


20 diagrams, cuts 


O'SMJSR 


0'AN,ISR 


3 diagrams, ISR 


<^SM,aUJSR 


'^AN,all,ISR 


20 diagrams, cuts, ISR 



Table 4: Cross sections and the corresponding physical assumptions under which they have 
been calculated. The subscripts SM, AN, on, off refer to Standard Model, non-standard 
TGC, on-shell and off-shell, respectively. 



and correspond to different physical asumptions. The flrst column refers to the SM and the 
second one to a non-standard TGC case (usually with only one of the CP-conserving couplings 
being different from its SM value). For the cross-sections labeled (Tcuts, cuts are applied mainly 
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to lepton and quark energies and angles in the laboratory frame: 



E,-,,,j>20GeV, |cos^,-,„,j| <0.9, cos^,_j<0.9, > 10 GeV . (19) 

The calculations were performed with input parameters as prescribed in the WW cross-section 
Working Group chapter. Results from the two programs agree within the MC errors. The 
particular case of daAN,off /dcosO (for the full phase space) has also been calculated by M. 
Bilenky in a semi-analytical method and full agreement with EXCALIBUR has been obtained for 
all CP conserving TGC. 
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Figure 3: Ratios of differential cross-sections at various levels of the simulation of the 4-fermion 

processes, (a) Ri = (TsM,off /<^SM,on, (b) R2 = (^SM,ISR/<^SM,off and (c) R3 = (TsM,all/<^SM,off,cuts- 



Different physical mechanisms could influence the angular distribution of the produced Ws 
and thus simulate the effect of non-standard TGC. Typical examples are shown in Fig. 
namely the effect of a finite W width, of ISR and of background graphs on da /d cos 6. ISR, for 
instance, lowers the available ^/s of the event and thus reduces the forward peak of the W~W^ 
production cross-section. In addition, the recoil of the W~W^ system against the emitted pho- 



ton further smears out the W angular distribution A similar effect, relative depletion of 



forward as compared to backward produced W~s can also arise from negative TGC parameters. 
This is evident from Fig. ^, where ratios of a non-standard da /d cos 9 and SM cross-sections 
are presented, both having been calculated under the same physical assumptions. Fig. ^(b) 
demonstrates the quantitative importance of this phenomenon. For final state electrons the 
background graphs, if not included in the analysis, could mimic a 6z of the order of —0.2. 
While the shape of the angular distribution da/dcosd for negative TGC parameters shows a 
trend similar to that induced by ISR, finite width or background graph effects, the normal- 
ization of the cross-section might provide some discriminating power, as do the decay angular 
distributions. Another very important message coming from Fig. ^ is that the sensitivity to 
the TGC remains the same at the different levels of the simulation (from on-shell Ws up to 
four-fermion production). Conversely, the influence of the various physics effects on production 
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and decay angular distributions is largely independent of whether or not non-standard TGC 
are present. 

We conclude that it is clearly important to account for and to correct the effects considered 
above in experimental analyses. We return to the effects of ISR and finite W width in Section 
5.2 where their neglect in TGC determination at LEP2 is quantified. In Section 6.2 we indicate 
how they contribute to the overall bias in a typical simulated TGC determination. 




Figure 4: Ratio of anomalous to SM differential cross-section, (a) <JAN,off /crsM,off (solid 
line), (Tan,isr/o-sm,isr (dotted line), (TAN,aii/o-sM,aii (dashed line), and aAN,aii,isR/ <ysM,aiiJSR 
(dash-dotted line) for = +0.1. (b) (TAN,off/crsM,off (sohd line), aAN,au/ crsM,aU muons 
(dash-dotted line) and electrons (dashed line) for aw = 0.2, 6z = 0.2, 6z = —0.2, aw = —0.2 
(bottom-top) and crsM,anl o'sM,off for muons (squares) and electrons (circles). 



4 Statistical techniques for TGC determination!^ 

Three different methods have thus far been proposed for the determination of TGCs at LEP2, 
— the density matrix method, the maximum likelihood method and the method of optimal 
observables. These methods are outlined in the following subsections and their application to 
common simulated datasets is compared. In devising these methods, two considerations have 
been borne in mind: first, — as will be elaborated in the next section — that it is advantageous 
to use as much of the available angular data for each WW event as possible; second, that the 

^^The experimental sections, ||-|[ have been coordinated by R. L. Sekulin 
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expected LEP2 data (a total of ~ 8000 events for an integrated luminosity of 500pb~ at 190 
GeV) will not be sufficient, for instance, to bin the data into the five angular variables appearing 
in the WW production and decay distribution (^) and subsequently to perform a fit- The 
studies reported in this section have been performed assuming that the final state momenta of 
the four partons from W~ and W^ decay have been successfully reconstructed from the data; 
the practical difficulties of doing this are discussed in section ^. 



4.1 Density matrix method 

In this method, TGC parameters are extracted from the data in a two-stage analysis. First, 
experimental density matrix elements and their statistical errors are determined from the an- 
gular distribution ( p!7D in bins of cos^; then the predictions of different theoretical models are 
fitted to the resulting distributions using a minimization method. The joint WW helicity 
density matrix elements Pxxx'x' are defined from (|l^) as the sums A, A)A^*(cr, A , A ) 

of bilinear products of production amplitudes and the dependence of the cross-section on the 
TGC parameters is fully contained in the complete density matrix thus evaluated. Similarly, 
by integrating over the observables of one W, single W density matrix elements Pxx' and Pxx' 
can be defined. 

The density matrix elements can be calculated in two ways: 

- Using the orthogonality properties of the W decay functions D^^> and D^^^ in (|^), density 
matrix elements can be extracted by integrating over the W decay angles with suitable 
projection operators. Thus, unnormalized density matrix elements of the leptonically 
decaying W in jjiu events can be found from the lepton spectrum as 

da(e+e- W+W') 1 f da(e+e- W+W- jjiu) ^ , , , , , , 

d^^;^ =B^J dcos^dcos^^.d^. ^.AO,<^d dcose.d^, 

(20) 

where Bwiu is the branching ratio for the jjiu channel, the angular variables are as defined 
in ([T3|) , ([I^) , with the decay angles and helicity indices now referring to the leptonically 
decaying W. Expressions for the normalized operators A^^^' are given in for example, 
Aoo = 2 — 5 cos^ 6'i projects out the longitudinal cross-section Poo^^^ of the leptonically 
decaying W. 



- In the second method 0, the production and decay angular distribution is expressed in 
terms of the density matrix elements and, in each bin of cos 6' , they are determined using 
a maximum likelihood fit to the distribution of the decay angles. 

Fig ^ shows some of the density matrix elements calculated from a sample of simulated events 
by the two methods as a function of cos 6' and fitted to the prediction of the Standard Model 
It can be seen that there is good agreement between the density matrix elements as calculated 
by the two methods, and with the fit to the Standard Model. 
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Figure 5: cos 6* dependence of density matrix elements pu and p_io for a sample of 2930 
simulated e+e~ — *• W~^W~ events at 190 GeV, calculated using the projection method (full 
circles) and the maximum likelihood method (triangles) and compared with the prediction of 
the Standard Model (fitted curve). 

4.2 Maximum likelihood method 

In this method, the distribution of some or all of the observed angular data is used directly in 
an unbinned maximum likelihood fit [0, in which parameters P, denoting one or more of the 
Lagrangian contributions (^, are varied to maximize the quantity 

In Cml = J2 P) - ^obs In jy p{fl, P)dn| , (21) 

where the sum is over events in the sample, Qi represents, for the i'th event, the angular 
information being used, p{fl,P) is derived from the cross-section (0), Nobs is the observed 
number of events, and the integral is over the whole of phase space. Many of the results 
shown here have been obtained using the method of extended maximum likelihood, in which 
the absolute prediction for the magnitude of the cross-section is also tested : 

IuCeml = ^lnp(a,P) -iV(P), (22) 

i 

where, for integrated luminosity L, the predicted number of events A^(P) in the sample is 
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L/^(f2,P)dO. 

It may be noted that, while in the evaluation of A^(P) in the absolute normalization 



of the cross-section must be used (as given in ([T7|) ), constant factors such as the flux factor 
may be omitted from the unnormalized expression J p{fl,P)dQ in (^). Furthermore, since 
for any event the probability p is proportional to the product of a phase space factor, which is 
independent of P, and a matrix element squared, | | , which contains the dependence on the 
TGC parameters, the sums over events in (|2lD and (^) may be replaced by In | Ai P), 
and the maximum of the likelihood function will be unchanged. While this replacement is trivial 
for the 2-body cross-section given by (p!2|), it is essential in the evaluation of the log-likelihood 
sum when the reaction is analyzed in terms of the 4-fermion processes, in which the phase space 
factor is different for every event. 

While the maximum likelihood method is able to use all the available angular information 
for each event, it has the disadvantage compared with a fit of being unable to provide a 
goodness of fit criterion. Nonetheless, the goodness of fit of a hypothesis represented by the 
likelihood function Ci{p) can be compared with that of C2{P) if the parameters p of Ci satisfy 
the condition p G P. Then the quantity —2ln{C^"'^/C^"'^), derived from the ratio of their 
likelihood functions, has a distribution |54|. This property has been applied to event samples 



generated with non-SM values of one TGC, Pi, and used to distinguish this hypothesis from a 
wrong one, when a different TGC, P2, is fitted to the data. — In general, a fit of P2 produces 
a result differing significantly from the SM value. Fig ^ shows the results of applying this test 
to the correct and wrong models in two alternative ways. In both cases, Ci is taken as the 
likelihood function when Pi varies; in the "same family" case (a), £2 is the likelihood function 
when both Pi and P2 vary, while, in (b), £2 describes a "composite" hypothesis, 

N 

C2 (Pi, P2; /5) = n If^P (^1) + (1 - /5) P (^2)] , (23) 

i=l 

where (5 is the probability that model 1, represented by the probability density function p{Pi), 
is correct, and Pi, P2 and (3 vary in the fit. It can be seen that a simple comparison between 
the values of these probabilities indicates the correct model for the majority of the cases. In 
addition, the absolute probabihty value indicates the goodness of the fit. 



4.3 Optimal observables method 

Optimal observables are quantities with maximal sensitivity |^ to the unknown coupling pa- 
rameters To construct them, a particular set of couplings Pj is chosen which are zero 
at Born level in the Standard Model (for instance, the TGCs defined by @)- Then, recalling 

This method has been used to search for CP vfolatfon m r+r^ productfon at LEPl, with a clear increase 
of sensitivity |55||. 
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Figure 6: Hypothesis testing using a) the "same family" and b) the "composite hypothesis" 
methods, for data sets of about 2500 jjiu events generated with TGC values deviating from 
the SM values by one to five times the expected LEP2 precisions. 



that the amplitudes for the four-fermion process are linear in the couplings, the differential 
cross-section may be written 

where Q represents the kinematic variables as before. Kinematic ambiguities, such as those 



described in table y, can readily be incorporated into (|2^) . The distributions of the functions 



are measured, and their mean values (Oi) evaluated|^. An example is shown in fig ^. To first 
order in the Pi, the mean values (Oi) are given by 



(O.) = (0.)o + ^Q,P, , (26) 
j 

from which the couplings Pj can be extracted because {Oi)^ and Cij are calculable given (|2^) 
and (^). From the distributions of the Oi the statistical errors on their mean values can be 

i^The functions ©^(Jl) for the TGC parameters used in ||] are available as a FORTRAN routine Q. 
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evaluated, the observables having been constructed to minimize the induced errors on the Pj. If 
the hnear expansion in the couphngs is good, the method has the same statistical sensitivity as 
a maximum likelihood fit. It can also be extended to incorporate total cross-section information 
in a manner analogous to the use of the extended maximum likelihood method discussed in the 
previous section. 



Optimal Observables - 190 GeV 




Figure 7: Distribution of Oq;^;^^^(cos^, (cos 6*/, 0/), (cos6'j, 0j)foidcd) for a large sample (50000) of 
simulated e'^e" — W^W~ events at 190 GeV. The experimentally determined mean value is 
to be compared with the expectation value of this observable in the SM, aw^j, = 0, used to 
generate the events. 



4.4 Comparison of methods 

In this section a comparison is presented of fits of the TGCs aw<t,, 0:34, and aw, defined in 



to common datasets generated with the PYTIIIA[57] Monte Carlo simulation program 



We precede this by mentioning the results of a comparison of the use of the maximum 
and extended maximum likelihood (ML and EML) methods, in which both of these methods 
were used in fits of the three TGCs to a large sample (50000) of events using first only the W 
production angle, and then the complete angular information (production and decay angles). 
The extra information contained in the EML method gave a substantial improvement (10%) 
in precision only in one case — the fit of a_B0, generally the least well determined parameter, 
to the production angular distribution. In the other fits the improvement was only ~ 1%. 
Similar conclusions have been obtained when applying the optimal observables method with 
and without total cross-section information. 

In the comparison of the density matrix (DM), EML and optimal observables (00) methods, 
the three analyses were applied to datasets at 175 and 190 GeV simulating both the expected 
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LEP2 statistics (~ 2000 events) and much larger statistics (50000 events). Sample results are 
given in fig ^, in which precisions obtained using the three methods in 1- and 2-parameter 
fits to the large dataset at 190 GeV are plotted. In all cases, the precisions obtained using 
the three methods are very similar when the same angular data is used in the fit. This can 
be seen in the figure, where the precisions from the EML and 00 methods, both of which 
used angular data cos^, (cos^;, (^i) and {cos 9j, (f)j)ioided, are almost identical. The DM results 
shown used the differential cross-section, ^j^^, density matrix elements Poo, Pi-i, Pio and p_]o 
of the leptonically decaying W, and the part symmetric in both polar decay angles of the 
transverse element ptt = Pii,ii + P-i-i,-i-i + Pii,-i-i + P-i-i,ii of the joint WW density 
matrix, representing somewhat less than the full 35 (CP-conserving) elements of the full joint 
density matrix. (Other density matrix elements can in principle be included in the analysis). 
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Figure 8: Comparison of TGC fits to a large sample of simulated events at 190 GeV using 
the density matrix (DM), maximum likelihood (EML) and optimal observables (00) methods, 
a): 1 s.d. precisions in 1-parameter fits to aw, (y-w<t> and as^- b): 95% confidence contours in 
2-parameter fits to {aw<t>i (^B<t>)- 

A difference between the EML or DM analyses and the 00 analysis can be seen in the 
2-parameter fit shown, where a second allowed region, remote from the SM region {awcp = 
0,0; _B<^ = 0) where the events were generated, is seen by the EML and DM methods. This 
effect is discussed in detail in ref. [^, where it is shown to arise naturally from the amplitude 
structure of WW production, and in particular from the fact that the helicity amplitudes are 
linear in the TGCs. It is not seen in the 00 results, because here the cross-section ( p4[) has 
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been linearized with respect to the TGCs about their SM values(^. 

In considering possible extensions to the analyses, two comments may be made. First, the 
EML and 00 methods could readily be used in a 4-fermion treatment by replacement of the 
matrix elements. The DM method does not lend itself to this adaptation, as the form ( ]T7| ) 
used in the projection of the density matrix elements assumes J = 1 for the two final state // 
pairs. Second, all three methods can in principle be adapted to the analysis of events with the 
experimental and other effects discussed later in this chapter; however, we have not made an 
assessment of the relative ease with which this can be done for the different methods. 

With the above points borne in mind, we can recommend all three methods for consideration 
in the analysis of LEP2 data. The studies reported in the following sections have, except where 
otherwise indicated, used ML or EML fits to obtain the results shown. 

5 Precision of TGC determination at LEP2: generator 
level studies 

In this section, the precisions to be expected in TGC determination from the anticipated LEP2 
integrated luminosity are summarized and an estimate of the biases and systematic errors 
accessible at generator level is given. 

5.1 TGC precisions in fits to simulated events 

Precisions in TGCs obtained from 1-parameter fits to simulated e'^e^ —>■ W~^W~ events at 176 
and 190 GeV are shown in table ^ and confidence limits in the planes of two of the three 
possible combinations of two of the parameters in eq. are shown in fig ^ Results are shown 
using various combinations of the angular data appropriate to each of the three final states 
jj^^y jjjj and iuiu, as indicated in table |^, as well as to the "ideal" case without angular 
ambiguities. For the first two channels (and for the "ideal" analysis), 1960 (2600) events were 
fitted at 176 (190) GeV; for the iuiu channel, 280 (370) events were used. These figures emulate 
the statistics anticipated from an integrated luminosity of 500pb~^ after experimental efficiency 
cuts of ~ 95%, 60% and 95% for the three channels respectively, and excluding leptonic decays 
into tUt-. The extended maximum likelihood method was used in the fits, and the events 
were generated and analyzed in the narrow W width approximation and without initial state 
radiation (ISR). In the analysis, the generated values of parton momenta were used, so that 
no account has been taken of the subsequent quark fragmentation nor of possible experimental 
effects. No kinematic cuts have been made on the data. The analysis reported here is therefore 
to be considered as an idealized one; the implications of the additional effects mentioned above 
are considered in detail in subsequent sections. 

^''An extension of the 00 method to incorporate second order terms in the parameters is under development. 
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Several conclusions may be drawn from inspection of the table and figure. As anticipated 
by the discussion of section substantial gains in precision are achievable by running at higher 
energy. Also, use of as much as possible of the available angular data serves to increase the 
precision and, in 2-parameter fits, to reduce the (quite pronounced) correlations between the 
fitted TGCs. The use of the jjjj channel, even with the angular ambiguities incurred by 
the inability to distinguish quark from antiquark jets, can be seen to provide a modest but 
worthwhile improvement in the overall precision attainable. Finally, the occurrence of a second 
region in the {aw(t>, dBct)) plane, remote from the Standard Model region (0,0) at which the 
events were generated but acceptable at the chosen significance level, has already been noted 
in the previous section. 
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Table 5: 1 s.d. errors in fits of (XB<i,-, OLw<i> and aw to various combinations of the angular data at 
176 and 190 GeV. The simulated data corresponds to integrated luminosity of 500pb~^. Details 
of the data samples are given in the text. 

In a first step towards a more realistic simulation of the data, some of the fits described above 
have been repeated using calculations corresponding to 4-fermion rather than WW production 
both in event generation and analysis. In so doing, contributions are included from the complete 
set of relevant diagrams and the finite W width effects ignored in the previous analysis are taken 



25 



e e 



1.25 - 



0.5 - 
-0.25 - 



-1 



-0.5 



S.2 

CD 

a 

1.25 
0.5 
-0.25 
-1 



COSl? 

COS-* 0, lej 




gO.8 















1 76 GeV 







sO.8 

0.4 


-0.4 
-0.8 , 



COSlJ' 

cos-tf e, lOjl 
cosiJ 0, 0, 




Figure 9: 95% confidence limits in the planes of 2-parameter TGC fits at 176 and 190 GeV, 
using various combinations of angular data, a), b), c), d): Fits to {aw<i>, o:b^)', e), f), g), h): 
Fits to {aw<i>, aw)- In the legend, the notation Qij implies a pair of decay angles {9ij, (pij) for 
W — > (leptons, jets) respectively, and | Qj \ implies the ambiguity cos^^ <-> — cos^j, 0j <-> (pj+TT 
incurred by the inability to distinguish quark from antiquark jets. In plots a), b), e), f), the 
angular data simulates channel jjiiy (and the "ideal" case, with no ambiguities); in c), d), g), 
h), it simulates channel jjjj. 
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into account. Using events generated with the ERATO |^6[ program corresponding to the 
expected statistics at 175 and 190 GeV, similar precisions to those shown above are obtained in 
fits of aw(j> and aw to angular data cos 6*, {cos 9i, (pi) and (cos6'j, 0j)foided|3- -'■'^ addition, in fits 
to a sample of Hlv events generated at 161 GeV corresponding to an integrated luminosity of 
lOOpb"^ (as suggested for the determination of the W mass from its threshold excitation [^), 
1 s.d. precisions of 0.18 and 0.43 were obtained in fits of awtj, and aw respectively. It is 
interesting to note that these values compare well with current experimental limits |2^, 26| , 
implying that TGC measurements from this exposure may also be of interest. This conclusion, 
however, remains to be tested when backgrounds and other experimental effects are included. 



5.2 Biases and systematic errors in TGC determination calculable 
at generator level 

It was pointed out in the previous section that the analyses presented there are idealized, in 
the sense that effects due to finite W width (unless a 4-fermion calculation is used), ISR, 
QCD and experimental reconstruction have been ignored. In this section, we consider the 
biases introduced in TGC determinations, first, if events generated with a realistic W mass 
distribution are nonetheless analyzed in the narrow width approximation, and, second, if ISR 
effects are also present, but ignored in the analysis. The discussion of the overall bias to be 
expected in TGC determination is pursued in the next section, where biases arising due to 
event selection and reconstruction are added to those discussed here. The systematic errors 
incurred both in the assessment of these biases and from other sources calculable at generator 
level are also estimated in this section. 

Figs. p!0|a) and b) show the effects of ignoring finite W width and ISR in the analysis of 
events generated with these effects included. Results are shown for several different generators, 
all operating in e"^e^ W^W~ (CC03) mode. It can be seen, first, that the bias incurred by 
neglect of ISR is greater than that from neglect of W width effects, second, that the biases are 
smaller when a fit involving more angular data is used, and, third (from b), that the biases are 
different for different values of a typical TGC parameter. Finally, we note that the overall bias 
is < the statistical error expected from LEP2 data. 

The systematic errors arising from these and other sources calculable at generator level are 
summarized, using a particular TGC fit as an example, in table BPl The first three entries 



come from the effects discussed above, the next two represent two different ways of expressing 
the uncertainty in the other electroweak parameters which are important in the evaluation of 



""^^A computational point may be made here: in the evaluation of the differential and total cross-sections 
needed in the likelihood expression (p2|), time may be saved by noting that, since the amplitudes for the process 
e+e~ /1/2/3/4 (or e+e" W^W~) are linear in the TGCs, an exact parametrization of the cross-section 
dependence on any one TGC may be found from a quadratic fit to its values for any three values of the TGC 
parameter. This procedure can be extended in an obvious way to fits of two or more parameters. 

^^The magnitude of some of these errors, in particular those arising from finite W width and ISR effects, 
depend on the angular data used in the fit, (c.f. fig p^ ). 
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Figure 10: Effect of ignoring finite W width and ISR in TGC fits, a): Results of fits of aw(i> 
to events generated with SM parameters at three energies using various generators. Left-hand 
plots: fit to cos6' only; Right-hand plots: fit to cos6', (cos^^i, 0i), (cos6'2, 02). b): as a), for 
EXCALIBUR events at 190 GeV, using cos6', (cos6'i, (pi), (cos6'2, 02), as a function of awij,- 
The legend for both plots is shown on b). 

the matrix element, and the final pair represent two independent uncertainties coming from 
machine and detector considerations. In any analysis which does not compare total cross- 
section predictions with the observed data, the second and last entries will not contribute to 
the overall uncertainty. It can be seen that, even when all the relevant entries are added in 
quadrature, the total is small compared with the statistical precision expected from LEP2 data, 
and we expect the larger component of the systematic error to come from uncertainties in the 
experimental effects considered in the next sections. 

In addition to the effects considered above, it is legitimate to ask whether colour recombi- 
nation effects among the two Ws could affect TGC measurements in the jjjj channel. It has 
recently been advocated that such effects may produce a shift of up to 400 MeV in My/ |59| . 
Therefore, by analogy with the effects of ISR, it may produce a bias in TGC measurements 
which would need to be accounted for, and, if not understood, would have an associated sys- 
tematic error. However, a preliminary study fS^] has indicated that the W production angular 



distribution, reconstructed from the hadronization products of generated jjij events, is little 



affected by application of the colour recombination models of ref |59|, and hence that it is un- 



likely that the shift in TGC values determined from the data in this channel will be significant 
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Source of uncertainty 


Uncertainty 


Systematic error 
in aw^ 


W width 


ATw = ±0.07 GeV 


±0.0004 


ISR 


Acrtot/(Ttot(e+e" W'^W~ + radiation) = ±1% 


±0.0013 


ISR parametrization 


Spread in Monte Carlo estimates 


±0.0020 


Mw 


AMw = ±0.18 GeV 


^0.0021 


sin^ 9w 


sin^ Ow = 0.226 (tree-level) sin' Ow = 0.231 


0.0029 


Beam energy 


Av^ = ±15 MeV 


^0.0002 


Absolute normalization 


±1% 


±0.0013 



Table 6: Systematic errors from various sources incurred in fits of aw(f> to angular data cos6', 
(cos 6*/, 0/), (cos6'j, 0j)foided at 190 GeV. The 1 s.d. statistical precision estimate for this fit 
from LEP2 data (c./. table |) is ±0.022. 



compared to the expected statistical error. 



6 Analysis of the jjev and jjjJ^i^ final states 



In the following we address some of the experimental aspects of the analysis of the e"'"e~ 
W~^W~ — >■ jjiu channel. In this section, we concentrate on the muon and electron channels, 
these being the cleanest and very similar in many respects. The tau channel is considered 
separately in the following section. For simplicity, the data are analyzed in terms of the five 
angles describing WW production and decay, by analogy with the generator-level analysis 
reported in section ^.1| . In its extension to a four-fermion treatment, also described in that 
section, the effect of the experimental selection and reconstruction procedures are expected to 
be the same. 



In section I^T^ we describe the efficiencies and purities obtained after the application of typical 
selection criteria and of kinematic constraints to the events. In the process of reconstructing 
and analyzing jjiu events, there are many experimental effects which can potentially bias the 
angular distributions, and hence the fitted values of TGC parameters. The scale of such effects 
is estimated in section |6]^, and in section we discuss briefly some methods proposed to 
allow for them in the analysis. The numbers presented result from a comparison of the work of 
several different groups and should be regarded as broadly typical of the four LEP experiments. 



6.1 Event selection, kinematic reconstruction and residual back- 
ground 



The jjiiy event selections used typically demand the following: 
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- that the event contains a minimum number, typically five or six, of charged track clusters; 

- that there is an identified electron or muon, or alternatively a high energy isolated track; 

- that the lepton has a momentum greater than its kinematic minimum, ~ 20 GeV; 

- that the lepton be isolated, by requiring low activity in a cone around the track (typically 
that the energy deposited in a cone of 100-200 mrad be less than 1-2 GeV). 

The effect of these cuts corresponds approximately to a fiducial cut in the centre-of-mass polar 
angle of the lepton of | cos 9iepton\ < 0.95. The acceptance for jets, which have some angular 
size, extends further but with falling efficiency. These numbers vary for specific detectors. 

The non-lepton system is then split into two (or more) jets using a conventional jet-finding 
algorithm. The following kinematic constraints can then be applied to impose energy and 
momentum conservation, and to improve the measurements using the fact that the system is 
overconstrained: 

ICfit: E^ = ^c™, Ep = 0, = 0; 

3C fit: In addition to IC, Mreconstructed = M]y for both W candidates; 

3C' fit: In addition to IC, Mreconstructed for both W candidates is constrained to a central 
value of MiY but is allowed to vary approximately within the W widtliQ. 

In the above, m,^ is the neutrino mass and M\y the W mass. A probability cut, typically 
of 0.1-1%, is applied to the constrained fit result. Typical efficiencies after these stages are 
shown in table ^ for centre-of-mass energies y/s = 175 and 192 GeV. The main loss is due 
to geometrical acceptance and lepton identification in the basic selection. The kinematic fits 
themselves are of the order of 90% efficient for such a probability cut. 

The background estimation was made using event samples, simulated with PYTHIA, of 
the final states WW (with neither of the bosons decaying to an electron or a muon), Z'j, ZZ 
and Zee. Also, contamination from 77 events, generated with TWOGAM were studied. 
Backgrounds from the last two channels were found to be negligible; those from the other final 
states are summarized in table |^. Contributions from the non-resonant graphs leading to the 
jjiu final state and containing TGCs have also been studied. It is found that, taken in isolation 
and ignoring interferences, they are rejected by the selection procedure. The main contribution 
to the WW background comes from events where one of the Ws decays into a tau and then 
into an electron or muon. Although this channel is sensitive to the TGCs, it will be seen in 
section |6.2| that the inclusion of such events into the analysis does not significantly bias the 
result. 

^^This is achieved by including either Gaussian approximations or true Breit-Wigner constraints in the fit 
procedure. 
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Table 7: Efficiencies and purities of the jjiiy sample at progressive stages of selection and 
kinematic fitting. 



Selection 


cos 9 


cos 6i 


Resolution 

(pl (rad) 


cos 9j 


(pj (rad) 


Before fit 


0.06-0.13 


0.11-0.17 


0.12-0.23 


0.13-0.19 


0.11-0.22 


After 3C' fit 


0.05-0.12 


0.07-0.13 


0.10-0.21 


0.10-0.17 


0.11-0.21 



Table 8: Resolutions on WW production and decay angles using simulated events at 192 GeV. 
The ranges indicate the spread of values obtained from different experimental simulations. 

Other approaches can be used instead of the selection procedure described above . In 
particular, if one wishes to avoid the use of the constrained fit, a cut requiring the missing 
momentum direction to be away from the beam pipe, typically cos^^ < 0.95, can be used to 
reduce the background from the ZZ and Z7 channels. In this case, an algorithm has to be 
applied to impose energy and momentum conservation. Nonetheless, in the rest of this section 
we adopt the 3C fits as representative of the efficiency and purity which can be achieved. 

The resolutions obtained for the WW production and decay angles before and after kine- 
matic fitting are shown in table |^. The values shown are averages over the whole fiducial 
region; however, in general, the resolutions depend upon the values of the kinematic variables 
themselves and, following kinematic fitting, they are correlated. It can be seen that a modest 
improvement in resolution is obtained, the main qualitative effect being due to the recovery of 
mis-measured events. 
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1-D 
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Statistical Precision 


±0.041 


±0.034 


±0.027 


±0.022 


(from table ^) 










Biases to Measurement: 










ISR and T\y 


-0.03 


-0.01 


-0.05 


-0.02 


Selection / Acceptance 


-0.06 


-0.02 


-0.03 


-0.03 


Reconstruction /resoln . 


-0.05 


-0.01 


-0.03 


-0.01 


Total 


-0.14 


-0.04 


-0.11 


-0.06 


Approximate additional 










bias due to backgrounds 










WW 


-0.005 


-0.002 


-0.003 


-0.002 


Zl 


+0.003 


±0.008 


-0.003 


±0.002 


zz 


-0.003 


-0.001 


-0.012 


-0.002 



Table 9: Biases in the measurement of aw4> estimated from studies of large samples of fully 
simulated events. In the last part of the table the additional biases due to residual backgrounds 
are shown. 

6.2 Systematic biases and statistical precision 

We now consider potential systematic biases, and the degradation of statistical precision due 
to experimental effects in the jjiiy channel. In this we include a) the neglect of ISR and r\y, b) 
experimental acceptance, c) reconstruction and detector resolution, and d) residual background 
contamination. The first item has been discussed in detail in section f).2\ , the result is included 
here for completeness. We use as example fits to aw<f> only. 

The overall bias due to a)-d) has been determined using a total of approximately 20,000 
simulated jjiu events at 175 GeV and 30,000 events at 192 GeV. A maximum likelihood or 
extended maximum likelihood fit was used, assuming in the analysis that the events originate 
from WW production with narrow W width and without initial state radiation. We emphasize 
that, since the purpose of this study is to show explicitly the scale of the biases, no corrections 
for the effects listed above have been applied in the analysis. 

The results are shown in table ^. The column labelled ID refers to fits using only the 
production angle cos^. The column labelled 5D refers to fits using the production and decay 
angles (with the angles of the hadronically decaying W folded to take account of the ambiguity 
described in table H). The bias due to ISR and Tw is derived as described earlier. The 
bias due to event selection and acceptance was determined by comparing fits to the generated 
angles before and after event selection, and the bias due to reconstruction and resolution was 
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determined by comparing fits to generated angles with fits to fully reconstructed angles. In the 
last part of the table the additional biases due to background are shown. However the reader 
should be aware that these were measured by adding small numbers of events to the sample, 
and in the absence of a systematic study should be considered to be very approximate. 

We conclude that the size of the biases from ISR and T^, acceptance and reconstruction 
are up to a few times the expected statistical error in the case of ID fits, and somewhat smaller 
when all the angular information is used. In order that these effects do not present a serious 
source of systematic error compared to the statistical error, they will eventually have to be 
understood and corrected for, incurring an error of less than ~ 10% of their values. 

Finally, we investigate the extent to which the statistical precision in TGC determination 
is degraded due to the effects mentioned above. The large simulated sample was divided 
into subsamples corresponding to the expected LEP2 statistics. The TGC parameter fit was 
performed on each sample, and the standard deviation of the spread of the results calculated. 
The precisions given for fits to generator level data for the jjiu channel in table ass ume an 
efficiency of 95%; thus the ideal precision in this channel is better by a factor VO.95 = 0.97. 
Taking this and the estimated experimental efficiency of 70% shown in table ^ into account, 
we expect a statistical degradation of ~ ±20% with respect to this ideal case. This is indeed 
observed, together with an additional degradation of ±10% to ±20% after application of the 
analysis procedure described above, showing the effect of the extra randomization from ISR, 
Tw and experimental effects. 

6.3 Strategies for allowing for systematic biases 

In the previous section the scale of the potential systematic bias due to detector and other 
effects was quantified. The simplest method of correction for such a bias is to determine its 
value for many simulated samples, subtract the mean bias from the experimentally measured 
TGC value and assign a systematic error on the basis of the width of the bias distribution 
and the experimental number of events. If the spread on the bias is large compared with 
the statistical error, this procedure will clearly be far from optimal. A second method is to 
use a Monte Carlo simulation to produce a correction function to map between "true" and 
"measured" values. This can easily be applied when fitting to a small number of variables, for 
instance to the cos 6' distribution alone, but is more difficult to apply in 5 dimensions simply 
because of the number of events required to characterize a 5D function in several bins per 
variable (unless corrections for each variable can be assumed to factorize). It has previously 
been shown at generator level that the precision is maximized by using all variables; it may 
however be that when systematic errors are taken into account the best overall precision is 
obtained by using a different strategy. 

It is nonetheless possible to formulate methods which take resolution effects into account 
in fits using all the kinematic variables. For instance, if the resolution/acceptance function 
for the variables ft is known, then the probability function p{Q, P) used for each event in the 
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maximum likelihood expressions (|2T| ) and (^) given in section |4.2| can be replaced by 



P) = / P(«.™.. P) X P(n,„„ ^ o_)dn,.„ (27) 

(where P represents the TGC parameters of the fit). The resolution/acceptance function p 
gives the probability that the true value O^tme would be reconstructed as 



"meas- 



There are several potential problems with the application of (P?]): (i) a 5-D integration 
is required; (ii) the resolution and acceptance functions will almost certainly not be simple, 
nor will they factorize; (iii) the correlations between angles must be known and included (in 



particular if kinematically fitted quantities are used). One suggested method uses fully 
simulated Monte Carlo events which are passed through the same events selection as data, in 
order to calculate the effective likelihood function. The variation of the TGC parameters is 
performed by reweighting the Monte Carlo events at their generated coordinates, while the 
comparison with data is performed at the reconstructed coordinates. This method can be 
applied for any fit dimension and can in principle take into account the effect of acceptance 
cuts, experimental resolution, any kinematic fitting procedure and background contamination 
in the data. 



7 Analysis of the jjru final state 

This channel requires special attention for two reasons. First, it comprises a sizeable part 
of the semileptonic WW decays and therefore could provide a useful addition to the available 
statistics and, second, it is a background mainly for the hadronic channel and therefore methods 
are required to reject it. 

In this study we consider only the hadronic decays of the r and describe criteria to select 
this final state. The resulting efficiency and purity expected for the sample and the resolution 
expected in the angular variables are presented. We find that an increase in the overall number 
of events selected for analysis in the jjiu channel of between 10 — 20% can be expected. 



7.1 Selection and reconstruction of jjri/ events 

To select jjru events, we make use of the characteristics of the r jet, namely small jet opening 
angle and low jet-charge multiplicity and of the global characteristics of the event, mainly 
missing energy and event acoplanarity. 

The signal for the jjru final state has a 3-jet topology, while the main sources of background 
(PFiy — > jjjj and WW Z'y{s) — >• qq'y) fall into the 4-jet and 2-jet topologies respectively. 
Thus the choice of the resolution parameter in a jet- clustering algorithm is quite significant. 
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Requiring at least 3 jets in the event, we find a r-reconstruction efficiency of 70 — 80% wliile 
only 30 — 40% of Z7 events survive. The clustering algorithm itself ensures isolation for the r 
jet. 

Jets from r decays can be distinguished from quark and gluon jets by the distribution of 
quantities such as the track multiplicity (total or charged), the maximum angle of any charged 
track in the jet to the jet axis, and the fractional energy of the jet contained within a cone of a 
specified angle (say, 0.1 rad) about the jet axis. A likelihood function based on such parameters 
has been constructed, giving a typical efficiency of about 70% with a rejection factor for quark 
and gluon jets close to 50. The charge of the r lepton can be estimated rather reliably from 
the total charge of the tracks in the jet (excluding those with momenta < 1 GeV/c from the 
sum in order to reduce the contribution from soft tracks from neighbouring jets). 

The T signal can be further enhanced by requiring that the event contains less than five 
jets and that the sum of the missing energy and the energy of the reconstructed r candidate 
should exceed \/s/2. This results in a selection efficiency for r events of about 90% with a 
rejection factor against the WW jjjj channel and against ZZ events of greater than 10. In 
addition, constraints on the polar angle of the missing momentum and the acoplanarity of the 
event can be imposed to reduce further the background from Z'j events. A rejection factor of 
10 is obtained while about 20% of the signal is lost. Finally, the very forward electromagnetic 
calorimetry can be used to detect ISR photon(s) in cases where they have not escaped in the 
beam pipe. 

It may be noted from the above that missing energy and missing momentum are key variables 
for the rejection of all types of background, and therefore the hermiticity of the detector is an 
important factor. 

The efficiencies and purities obtained for jjru events from a sample of simulated events at 



192 GeV are shown in table |I0| The background from the jjii' channel stems mainly from 
inefficiencies in muon detection in the simulation used, and some improvement may be possible 
here. The application of a 2-constraint kinematic ffi0 can also be seen to provide background 
rejection, with a small decrease in the r selection efficiency. 



Selection 


Efficiency % 


Background % 








WW - 


jjjj 


WW jjiu 


ZZ 


total 


No fit 


35 - 45 


4 - 6 


4 - 


8 


5 - 8 


- 2 


13 - 24 


2C ffi 


32 - 42 


- 2 


2 - 


5 


5 - 8 


- 


7- 15 



Table 10: Typical efficiencies and purities for the jjriy channel with no kinematic ffi and with 
a 2-constraint kinematic fit. 

^^The 2C fit imposes energy and momentum conservation and constrains the jj and ti^t systems to have the 
W mass, leaving the momentum of the neutrino from W decay and the t energy as free variables (with a lower 
limit on Er given by the visible energy of the r decay products). 
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An improvement to the kinematic fit may result by constraining the r momentum, using the 
fact that the direction of the r can be accurately estimated from the combined momentum 
of its visible decay products, so that the r energy can then be computed from the W decay 
kinematics p3 . 



7.2 Resolution in reconstructed quantities 

The resolution in the centre-of-mass polar and azimuthal angular variables of the r, evaluated 
using 2-Gaussian fits to the differences between reconstructed and generated values, is of the 
order of 5 mrad in 75% of the events, and is not changed much by the kinematic fit. The 
energy of the original r can only be estimated at a level of AE/E = 0.15 with no kinematic 
fit, but after the fit has a resolution AE/E = 0.05 in 80% of events. The resolutions in the W 
production and decay angles, evaluated after the kinematic fit, are found to be Acosd = 0.11, 
AcosOj. = 0.13 and A^,- = 250 mrad respectively. 



7.3 TGC determination from jjriy events 

The precision with which TGCs can be determined from jjru events has been investigated 
using a sample of 937 fully simulated events, generated at 192 GeV with finite W width and 
ISR, corresponding to an integrated luminosity of 500pb~^. Of these events, 390 survived the 
selection and reconstruction procedures described above. The parameter aiy</, was fitted to 
the cross-section (|17D {i.e. in the narrow width, no ISR approximation), using the extended 
maximum likelihood method described in section and folding over the 2 ambiguous solutions. 



The 1 s.d. precision in awcf> was found to be ±0.06 with estimated biases of —0.04 from the 
neglect of W width and ISR, —0.025 from the effects of reconstruction, and +0.03 from the 
presence of background events. 



8 Analysis of the jjjj final state 

The advantage of the high branching fraction of this channel is somewhat reduced by exper- 
imental difficulties associated with the purely hadronic nature of the final state. Background 
rejection in the four-jet channel is difficult, since no high-energy charged lepton is present to 
tag one W as in the semileptonic case. The largest background is expected from the high 
cross-section channel e"'"e~ — >• Z/7* — qq(7) leading to multi-jet final states. Also, since the 
decay modes of the two Ws are both hadronic, the problem arises of selecting the correct pairs 
of jets to form the two Ws and of assigning their charges. 

In the following we suggest an analysis of the jjjj channel, including event selection, jet 
reconstruction and kinematic fitting, and indicate the expected efficiency and background levels. 
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A section is devoted to jet- and W-charge tagging. We then discuss the determination of TGCs 
from the selected events, and draw conclusions on the sensitivity of the four-jet channel. 



8.1 Selection of events and reconstruction of 4 jets in the final state 



The general criteria for the selection of jjjj events are based on the fact that the hadronization 
of four quarks gives rise to a high multiplicity of particles in the final state, and to a large visible 
energy. Other types of events with hadrons in the final state can have similar characteristics, 
mainly the jjiu channel and the reactions e~^e~ qq'y with Mqq >120 GeV0 and e~^e^ 
ZZ qqq'q' ■ The first two reactions can mimic 4 jets when gluon radiation has occurred. 

The following variables were typically used to select jjjj events: 



- A large multiplicity of particles in the detector {N charged > 25, or N charged + ^neutr > 
25 — 40). This cut helps to reject jjiu and QCD background, where the observed hadrons 
originate from a smaller number of initial quarks; 

- Small thrust and/or large sphericity (T < 0.9 — 0.97 or S > 0.1). The QCD background 
generally consists of two back-to-back jets (T ^ 1, S* — 0), while the WW hadronic 
decays are more isotropic. However, the discriminating power of these variables becomes 
smaller as y/s increases; 

- Large total visible energy (charged -|- neutral); 

- Small missing energy {Emiss < 40 — 50 GeV). Large missing energy and momentum are 
associated with the neutrino in leptonic W decays, and with an undetected high energy 
photon in qq'y events. 



Events from the jjii^ channel can also be suppressed by requiring that no energetic isolated 
track or high energy identified lepton be present. The efficiency of the selection criteria at this 
stage is typically around 80% and the purity of the surviving sample is around 60%. 

After the cuts described above, the final state particles are grouped into four jets. For this 
purpose, several clustering algorithms have been tried, which fall into two categories, namely 

DURHAM 13 or 



transverse momentum-based clustering, such as LUCLUS pSl, PUJET4 p6 



GENEVA [^, and scaled invariant mass squared clustering, such as JADE |6^. Comparative 
studies have shown that differences are contained to within about 3%, with the algorithms based 
on transverse momentum reproducing the initial parton directions somewhat better, leading to 
better jet definition and hence better resolution in invariant mass. 



^^Events with a lower invariant mass of the qq system correspond to radiative return to the Z*^ peak and can 
be easily rejected either because the photon radiated from the initial state is detected or because the missing 
momentum associated with it is very high. 
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Further rejection of background can be achieved by apphcation of the following additional 
cuts to the reconstructed jets, leading to a jjjj purity of ~ 80%: 

- Minimum number of particles inside each jet (2 to 5); 

- Minimum angular separation between jets (20") ; 

- Minimum energy of the 2 least energetic jets (15-20 GeV); 

- Minimum jet-jet invariant mass ( Ycut = 0.002^/s). 



8.2 Kinematic fitting 

The kinematic fit is used as a tool to improve the resolution on measured quantities by impos- 
ing external constraints. For the jjjj channel, the measured quantities are the energies and 
polar and azimuthal angles of the four reconstructed jets (and, for massive jets, their invariant 
masses). The external constraints which can be imposed are as follows: 

1) energy-momentum conservation (4C), 

2) as 1), plus equality of the two reconstructed invariant jet-jet masses (5C), or 

3) as 1), plus equality of the two reconstructed invariant jet-jet masses with Mw (6C). 

The importance and the limits of kinematic fitting have been discussed in previous sections 



of this report, and technical details can be found in references |66|, [70[|. As in their application to 
TGC determination in the jjii' channel (see section ^ the second and third constraints can be 
imposed without fear of introducing biases, as they would if applied to W mass determination. 
Nonetheless, a comparison of results using different constraints has shown that there is negligible 
gain in going from the 4C fit to the 5C or 6C fits, and the results given below have used a 4C 
fit, followed by cuts on the invariant masses of the jets assigned to each W. In order to choose 
the best pairing of the four jets into two Ws , kinematic fits are made to each of the three 
pairings, and that with the largest probability is taken as the correct combination. 



8.3 Results in efficiency and resolution 

After additional cuts on the fitted quantities to reduce background contamination, the effi- 
ciencies, purities and remaining background content of selected event samples generated with 
different detector simulations and at two centre-of-mass energies are as shown in table p!l| . 



true 
et 



The resolutions in the radial and azimuthal jet angles 6 jet-, 4>jet and the resolution AEjet/Ej 
in the jet energy can be estimated by comparing each reconstructed jet with the closest gen- 
erated quark direction. They show little dependence on the centre-of-mass energy and on the 
different detector simulations. Results for the resolutions in jet energy and in the reconstructed 



W production angle cos^ for simulations at 192 GeV are shown in table ^ It can be seen 
that the kinematic fit substantially improves the resolutions in the variables shown (by a factor 
of between 30 and 50%). However, it has less effect on the jet angular resolutions, which are 
typically between 20 and 30 mrad for about 2/3 of the selected events. 
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y/s = 175 GeV 


v/s = 192 GeV 


Efficiency (%) 


54 - 59 


52 


Purity (%) 


92 


90 


Background(%) 






e+e^ — i> qqj 


8 


8 


e~^e~ — >■ ZZ qqq'q' 





2 


e+e" ^ WW jjiu 









Table 11: Efficiency and purity of samples of events selected with the cuts described in the text 
at two centre-of-mass energies. 





Before kinematic ffi 


After kinematic ffi 




0.12 


0.08 


Acos^ (mrad) 


50.0 


40.0 



Table 12: Resolutions in jet energy and W production angle before and after the kinematic ffi 
at 192 GeV. 



8.4 W charge assignment 



The ambiguities in angular data arising from the inability to distinguish quark from antiquark 
jets in W decay have been hsted in table ^, and the generator level studies simulating the 



jjjj channel described in section ^TT| were made using distributions folded in both production 
and decay angles. In order to attempt to resolve the ambiguity on the production angle, a jet 
charge can be defined by weighting the charge Qi of each particle assigned to the jet with some 
function of its momentum, 



Different weight functions have been tried, based on transverse momentum, on rapidity, and on 



a power of the momentum |171| , |72| , ffSj |f^ . It appears very difficult to identify the charges of each 
individual jet. But, since the separation between the charges of the two Ws is equal to 2, one 
can more easily distinguish the W~ from the W^ and therefore determine the production angle 
in the lab frame. The charges of the two jets assigned to one W on the basis of the kinematic 
fit are therefore added together to evaluate the charge of the W . The fraction of selected events 
where the charge is correctly assigned is found to be 80%. No significant difference among the 
various weight functions was found. The W charge identification implies a gain in sensitivity 
in TGC determinations. 
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8.5 TGC determination from jjjj events 



The precision obtained in TGC determination after application of the procedures outhned 
above has been estimated using a fully simulated sample of 2292 events at y/s = 192 GeV, 
corresponding to an integrated luminosity of 500 pb~^. Two types of fit were performed to 
the observed angular distributions, namely, a fit to the production angle cos^ only, and 



an unbinned maximum likelihood fit 
folded W decay angles (cos6' 



fas described in section 



to the production angle and 
In both fits, the ambiguity in 
production angle was resolved using the jet charge assignment. Precisions obtained in fits to 



n 



) folded, 



fcos6' 



j2 1 r'j2 ) folded • 



the TGC parameters aw^p and are shown in table |13 



Fitted parameter 
value 


Fitting method 


method 


Maximum likelihood method 




0.04 


0.02 


aw 


0.07 


0.04 



Table 13: 1 s.d. errors in fitted values of parameters awct> and aw to a sample of 2292 fully 
simulated jjjj events at 192 GeV. fits were made to the production angle only and maximum 
likelihood fits to production and folded decay angles. 

The data used in the fits were generated according to the Standard Model using PYTHIA, 
with Tw = 2.1 GeV and with ISR. The theoretical expectations |^ were calculated with 
= and without ISR. In these conditions, a biased result is expected, as indicated from 
the results shown in fig ^ (section |5.2| ). In addition, experimental biases due to the selection 
and reconstruction procedures are to be expected, as found in the analysis of the jjiiy channel 
and discussed in section In the case of the jjjj channel, the angular distributions are quite 
severely distorted by bad association of pairs of jets to the parent W and by wrong W charge 
assignment, and the resulting biases can easily simulate an anomalous TGC. The results shown 
for the fit to the production angle only include the effect of the application of a procedure to 
correct for the bias. Although based at present on the use of very limited Monte Carlo statistics, 
the fitted central values are found to remain within ~ la of the SM values after application 
of the correction. However, a full study of the biases in this channel and the development of 
methods to correct for them in fits using several angular variables have yet to be carried out. 



9 Analysis of the iuiu final state 

The analysis of the channel in which both Ws decay leptonically presents particular problems. 
It is the least statistically significant final state (with branching ratio ~ 11% for / = e, /x, r), the 
missing neutrino momenta imply that the W direction cannot be determined unambiguously, 
and, if one or both of the Ws decay into rz7^ or its charge conjugate, the presence of the 
extra neutrino from r decay makes it impossible to reconstruct the event, reducing the useful 
branching ratio of such events to around 5%. On the other hand, the knowledge of the W 
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Cut 


]\j 9 

' leptons ^ 

leptons G {e, /i} 


P:p''' > 1.5 GeV 


M,y<Mz-Tz, 

' ■ Mfp> Mz + Tz 


Recon- 
struction 


(GeV) 


175 


190 


175 


190 


175 


190 


175 


190 


Efficiency (%) 


82.7 


80.1 


82.6 


79.9 


79.2 


77.8 


58.7 


58.5 


Purity (%) 


9.70 


9.35 


25.2 


24.5 


31.9 


30.9 


88.4 


80.8 


Background (%) 




















86.0 


87.0 


63.6 


65.9 


53.4 


57.7 


6.33 


14.1 


ZZ 


0.4 


0.5 


1.22 


1.32 


1.60 


0.89 


0.32 


0.28 


WW ^ Ivrv 


3.83 


3.16 


9.99 


8.28 


13.1 


10.5 


4.93 


4.75 



Table 14: Efficiencies and purities in selection of £i/£z/ events at 175 and 190 GeV. 

charge and the small reconstruction errors of leptons favour this channel in contrast to the 
4 jet channel. In this section the usefulness of the purely leptonic decay channel for TGC 
determination is discussed in the light of these considerations. 



9.1 Selection of events 

The event signature is very simple: two leptons and large missing energy. This makes the 
channel easy to identify, but the background contributions, chiefly from Z7, are high. Also, iuiv 
events with one or two leptonic r decays {BRr-,e,^l ~ 35%) constitute a possible background 
of about 1.8% of the total number of WW events. The typical selection criteria used for iuiu 
events aim at reducing these backgrounds by requiring large missing transverse momentum 
and, for equal flavours, that the mass of the lepton-lepton system should not be close to Mz- 
In addition it is also necessary that physical solutions to the reconstructed neutrino directions 
exist - this turns out to give the strongest background rejection. 

For purely leptonic WW events the momenta of the 2 neutrinos are unknown. However, 
in the absence of ISR and for fixed Mw, we have six constraints allowing the momenta of the 
neutrinos to be reconstructed 0. The quadratic nature of these constraints results in a two- 
fold ambiguity, corresponding to flipping both neutrinos with respect to the lepton-antilepton 
plane, hence only affecting cos 6', (pi, and 02, while leaving cos^i and cos6'2 unchanged. 

The efficiencies and purities after each stage in event selection and reconstruction are shown 
in table for fully simulated events generated with ISR and finite width. It can be seen that 
the required existence of solutions to the six constraints provides a very strong background 
suppression. However, it is important to note that the solution of these equations requires the 
use of all the kinematic information available in the event, leaving no possibility, for instance, 
of accounting for ISR or finite W width effects. Thus, with these effects included, no solution 
is found at generator level for about 20% of the events. 
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9.2 TGC measurements from ^v^v events 



The precision with which TGCs can be determined from events has been investigated 
using samples of fully simulated events, generated at 175 and 190 GeV with finite W width 
and ISR, corresponding to an integrated luminosity of 500pb~^ The parameter aw,t> fitted 
to the cross-section (pTf) {i.e. in the narrow width, no ISR approximation), using the extended 



maximum likelihood method described in section and folding over the 2 ambiguous solutions. 
The 1 s.d. precision in aw<i> was found to be ±0.15 at 175 GeV, with estimated biasesQ of 
—0.04 from the neglect of W width and ISR, —0.05 from the same sources plus the effects of 
reconstruction, and a combined bias of —0.07 when, in addition, background events are added. 
At 190 GeV the precision in aw(t> was found to be ±0.09 and the same biases —0.04, —0.13 and 
—0.21, respectively. 

Taking into account the small number of events (~ 220) in the sample, it is clear 
that the sensitivity to TGCs is highly preserved in this channel, despite the two-fold ambigu- 
ity. However, it is clear that, due to the very limited statistics, they will have to be used in 
combination with events from other decay channels. 



10 Other Anomalous Couplings and Other Channels 

10.1 Constraints on WWj Coupling from e^e~ Vu^ 

The W^W^ production process suffers from the drawback that both W^W'^ and W^W~Z 
couplings contribute and it is not easy to disentangle the various contributions. However, 
there does exist a process, namely e"^e^ i^^l-, which allows us to concentrate solely on the 
W^W~'~^ vertex. The matrix-element for Vvy production in terms of the WW'j TGCs k^, 
in (|l]) has been calculated in Ref.[|75[. In the numerical analysis we set acceptance cuts of a 



minimum photon angle of 20° and transverse momentum of 10 GeV. To increase the sensitivity 
to anomalous couplings the background from the Z exchange graphs, e^e~ — >■ Z'y — >■ Vwy, 
is eliminated by requiring the energy of the photon to be at least bTz away from the energy 
corresponding to the Z'j final state which essentially amounts to an upper limit on photon 
energy of 53 GeV. With these cuts the cross- sect ion0 for the standard model is 1 pb at ^/s = 175 
GeV, which still leads to an appreciable number of events at design luminosity of 500 pb~^. 
Cross-sections for non-standard TGC, within these cuts, differ by less than 0.1 pb for |Afi:| 
and/or |A| < 2, so not much sensitivity is expected from the total cross-sections alone. Looking, 
however, at the deviations of the differential cross-sections from the standard model predictions 
one can set some limits on the parameters. We consider a fit to SM data, adding in 

^°Due to limited statistics the statistical errors on the results from which the biases are estimated are of the 
same order as the biases themselves, but since the samples are correlated the statistical error of the biases are 
expected to be smaller. 

^^We have not included effects of initial state radiation. 
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quadrature a relative systematic error of e = 0.02. In Fig. |ri| we show the contour plots for the 
distributions as functions of Ak^ and as extracted from a) the energy, b) the transverse 
momentum distributions of the photonQ. We used equal size binning with 17 and 16 bins for 
the two cases respectively. This process is, in general, more sensitive to Ak^ than to A^. It 
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Figure 11: contours in the A/t^-A-y plane as derived from (a) energy and {b) transverse 
momentum distributions, for ^/s = 175 GeV and integrated luminosity of 500 pb~^. 



is thus somewhat complementary to Tevatron bounds which are derived from W'j production. 
While quantitative improvements on the constraints may be made by considering two-variable 
distributions or by adopting maximum likelihood methods, these would still not be competitive 
with those deduced from W^W~ production. However, the Vv^ channel isolates the WW'y 
couplings and probes them in a different q"^ region. Therefore it complements the information 
obtained from l^-pair production. 



10.2 Anomalous couplings Q 

While the measurement of WW'j and WW Z couplings at LEP2 has deservedly received con- 
siderable attention, it is also important to search for couplings between the neutral gauge 
bosons |]75|, [TTI] . For the trilinear ZV^7 vertex {V = 2', 7) the most general vertex function in- 
variant under Lorentz and electromagnetic gauge transformations can be described in terms 
of four independent Q dimensionless form factors]?^, denoted by hY , 1=1,2,3,4. The parts 
of the vertex function proportional to and are CP-violating while those involving the 
other pair of form factors are CP-conserving. As is well known, all Z7 form factors are zero at 
the tree level in the SM. At the one-loop level, hY and hY are zero while the CP-conserving 

^^The angular distributions are less sensitive to the anomalous couplings. 

We are grateful to Ulrich Baur for making his Zj event generator available to us. 
24 As for the WWV TGCs of Eq. (|l|), constraints on the different hY can be obtained from restriction to the 
lowest terms of a gauge-invariant expansion in l/A^p. 
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form factors are nonzero but too small to lead to observable effects at any present or planned 
experiment. Observation of Z7 couplings would, therefore, signal physics beyond the SM. 

We have carried out a generator-level study to estimate the sensitivity at LEP2 to anomalous 
couplings. Following reference |]7B[, the form factors were parameterized as /i]^ = + 
(P^/Ay))"'» where P is the effective center-of-mass energy, and hjj^. Ay, and are parameters. 
For comparison with present limits on Z7 couplings, we chose ni = 71,3 = rig = 3.0 and 
n2 = n4 = rig = 4.0. Two channels, e~^e~ —>■ fi^fi~'~f {fifj^'y) and e~^e~ —>■ uu'j (I7), have been 
studied in detail. At LEP2 energies it turns out that the I7 channel is much more sensitive 
to anomalous Z7 couplings than the fifi'j channel. This is due mainly to anomalous couplings 
being dominated by the case where the detected photon recoils against a resonant Z and that 
r{Z —>■ vv) = Qr{Z fi^fi^). Below we thus report on the sensitivity expected from the I7 
channel alone. 

Experimentally, anomalous couplings in the I7 channel would populate the same energy 
range as "radiative return" to the Z pole through initial-state radiation (ISR), namely, the 
reflection of the Z pole centered on Eq = {s — m|)/(2-v/i). Unlike photons from ISR, however, 
photons from anomalous couplings are distributed almost uniformly in solid angle. In our sen- 
sitivity analysis, which employed event counting rather than fits to distributions, we therefore 
required (a) the photon energy to be within 10 GeV of Eq and (b) \cos9y\< 0.8 in order to 
maintain good acceptance for anomalous couplings while suppressing the background from ISR. 
For I7 events passing these cuts, a combined trigger and reconstruction efficiency of 90% was 
assumed. 

Figure |T3(a) shows the ZZj couplings that would be excluded at the 95% C.L. for y^=175 
GeV and 500 pb~^ assuming that the SM expectation is observed^. The limits are shown for 
two different values of Az to provide some indication of how much they depend on the particular 
choice of parameter values. Limits on these couplings have been published recently by L3||79| 



CDF[80|, and D0|8l|]; the L3 and CDF limits are also plotted. It is evident that the expected 



sensitivity of LEP2 is comparable to the combined sensitivity of searches by LEPl and Tevatron 
experiments. Figure [l2|(b) shows the corresponding estimated sensitivity to anomalous Z77 
couplings. As can be seen from comparison with the limits from CDF||8y] (competitive limits 



have also been published by D0[^), LEP2 is expected to extend considerably the sensitivity 
to Z77 couplings. 

The sensitivity to anomalous Z'j couplings increases rapidly with center-of-mass energy, 
the effect being more pronounced for sensitivity to and , which correspond to dimension-8 
operators compared to dimension-6 operators in the case of and . For example, sensitivity 
to hjf) (hlo) is improved by about 25% at 192 GeV, even with a smaller integrated luminosity of 
300 pb~^. Although backgrounds are expected to be more severe, analysis of the event sample 
consisting of hadrons and an isolated, energetic photon may provide another way of significantly 
increasing sensitivity to Z7 couplings. 

The effects of QED corrections on LEP2 sensitivities are not reflected in Fig. |2[ These corrections reduce 
the sensitivity to anomalous Z7 couplings but by less than 10%. 
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Figure 12: Estimated LEP2 sensitivity limits (95% C.L.) to anomalous Z7 couplings and 95% 
CL. limits from present experiments. The LEP2 estimate is for -^5=175 GeV and 500 pb^^. 
See text for explanation of the parameters. 



10.3 Constraints on gauge boson interactions from e ^ qq ,11 

The description of NP for e~e~'" qq ,11 in terms of dimension 6, purely bosonic, SU (2) x U{1) 
gauge invariant operators necessitates the consideration of the interactions 

^NP = -^^{D,W^,)-{D^W') - ^-^{d,B^,){d^B^n 

NP NP 

- <^^B,,f-W^<^ + ifL(D,,$)t$$t(D^$) , (29) 

in addition to the ones mentioned in Section |2.1.1| . Such interactions affect the gauge boson 
propagator at the tree level and are thus rather strongly constrained by LEPl measurements. 
Nevertheless LEP2 can significantly improve these constraints, particularly for the first two 
terms in ([25|) which give a q^ contribution to the gauge boson propagator It has been 



remarked in that if the physical quantities measurable in e~e+ QQ , at LEP2 are 
expressed in terms of Z-peak observables, then the aforementioned contribution allows the 
remaining anomalous dependence of the results to be described in terms of only the two cou- 
plings fow and foB- Thus by looking at cr-^Q^^YOYis ' ^^^^ strong constraints on these 
couplings should be possible; (see Fig. |13|). 
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Figure 13: Sensitivity of LEP2 to fnw and fcB from e"e+ qq ,11 at y/s = 175 GeV and 
500 pb^^ (one experiment). Constraints from (Thadrons (solid lines); cr^+r (dashed lines); A'p'^ 
(dash-dotted lines); ah (dotted lines); global fit (solid ellipse). A^rp = 1 TeV is assumed. 



10.4 Higgs anomalous couplings 

Anomalous couplings could also arise for the Higgs interactions with itself and the gauge bosons. 
In fact, dynamical considerations indicate that it is easier to generate anomalous couplings for 



the Higgs rather than for the gauge bosons [|1^, [T3[. The dimension 6, SU{2) x f/(l) invariant, 
CP conserving interaction is 



Cmp = ^ ($^$ --){dW ■ W,, + ds B'^-' B^,) + 9^($t$)a^($t$) . (30) 



The first two of the above terms generate Higgs-gauge boson interactions while the last one 
induces anomalous Higgs interactions through a renormalization of the Higgs field. 

As in section |2.1.1| , unitarity can be used to associate to any given value of each of these 
anomalous couplings the largest allowed scale A^ where New Physics generates it. For the first 
two operators these relations are 

104.5 104.5 (^^"^ 



d \ \ for d > , d- ) ( for < , (31) 



l + 6.5(^ l-A'"^- 



w 
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195.8 (f^ 

1 +200 [^y 



for (Ib > , d 



B 



195.8 (fg^ 
1 + 50 [M^] 



for c/r < . 



(32) 



Thus, for A[/ = 1 TeV, the largest allowed values are d ~ 0.4 or —1 and ds — 0.6 or —1. 

The above anomalous Higgs couplings may be studied at LEP2 through the processes 
e~e^ ZH, provided mn < \/s — Mz, or via e~e'^ 'yH if rriH < y/s. Considering tree 
level anomalous contributions and restricting to cases where only one of the operators above is 
active [p3| , ^ p5| , |86| , we get the results given in the figures below. Thus, from Fig. [l^a, pre- 
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Figure 14: Distribution of Higgs production angle, da/dcos9, for (a) HZ production at uih 
80 GeV and (b) H'y production at rriH = 120 GeV. 



senting e e+ — ^ ZH, we deduce observability limits |/<^2| — 0.01 and \d\ 
corresponding to Au ~ 6 — 7 TeV ( Au ~ 5 TeV) for uih — 80 GeV. 



0.015 (I^bI ~ 0.05) 



More striking is the process e~e^ 'yH which is unobservable at LEP2 in the SM |55| , 
but may become observable in the presence of NP interactions for the Higgs. A sensitivity to 
\d\ ~ 0.05 or Idsl — 0.025 should be possible from this process for tjih ~ 80 GeV, which means 
testing NP scales up to 3 and 7 TeV, respectively[|86[ . 



11 Conclusions 



Experiments at LEP2 will allow a precise direct measurement of the most immediate conse- 
quence of the non-Abelian character of the electroweak bosons, the TGC of the W to the 
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photon and the Z. Various channels can provide information on non-standard interactions in 
the bosonic sector. The process e^e^ ff determines obhque parameters which are comple- 
mentary to LEPl results. Z7, HZ and Hj production allow one to search for non-standard 
boson couplings in the neutral sector. e~e"'" i/Pj is marginally sensitive to the WWj cou- 
pling in isolation. However, the most important process is clearly e~e'^ W~W~^ or its 
generalization, 4 fermion production. 

Of the various decay channels, the semileptonic modes W~W^ jjG^-, jjf^i^ will provide 
the most precise individual measurements of TGCs, since high statistics and almost complete 
information on the decay distributions are combined. Of particular importance is the iden- 
tification of the W charge which is needed to measure the full production angle distribution 
da/dcosd. Also, the decay angular distributions and their correlations with each other and 
with the W production angle are needed to resolve the correlations between different TGCs to 
a maximal extent. 

A priori, the jjjj final state provides incomplete information on the W charges. How- 
ever, correct charge assignments at the 80% level can be obtained by determining weighted jet 
charges, providing potentially valuable additional information in TGC determination. While 
more limited in statistics, the leptonic channel, is particularly clean, and the jjrp channel 
will also be of use in TGC analyses. 

Using jjev and jjiiu data alone, measurements of particular TGC parameters at ^/s ~ 
192 GeV appear possible at generator level with a precision of ^ ±0.02 for an integrated 
luminosity of 500 pb~^. The effects of ISR and finite M^- width and the apphcation of exper- 
imental selection, acceptance and reconstruction procedures lead to a degradation estimated 
at ^ 30 — 40% in the precision, and to a systematic shift which is a factor 3 larger than the 
statistical error, but our studies indicate that this bias can be corrected. For more general 
TGCs, considerable cancellation between different parameters is possible, resulting in weaker 
bounds. It is for this case that information from the full five-fold angular distribution of 
production and decay angles or its generalization to 4-fermion final states becomes particularly 
important. 
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